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Abstract 

We review the apphcation of non abelian discrete groups to tlie theory of neu- 
trino masses and mixing, which is strongly suggested by the agreement of the Tri- 
Bimaximal (TB) mixing pattern with experiment. After summarizing the motiva- 
tion and the formalism, we discuss specific models, based on A4, 6*4 and other finite 
groups, and their phenomenological implications, including lepton fiavor violating 
processes, leptogenesis and the extension to quarks. In alternative to TB mixing 
the application of discrete fiavor symmetries to quark-lepton complementarity and 
Bimaximal Mixing (BM) is also considered. 
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1 Introduction 



Experiments on neutrino oscillations, which measure differences of squared masses and 
mixing angles P-Ej have established that neutrinos have a mass. We refer in particular 
to ref. [T] for an introduction to the subject, the main results, the basic formalism and all 
definitions and notations. Two distinct oscillation frequencies have been first measured 
in solar [3^9] and atmospheric [TOl[Tl] neutrino oscillations and later confirmed by exper- 
iments on earth, like KamLAND [12], K2K [13], MINOS [H] and OPERA [15]. A signal 
corresponding to a third mass difference was claimed by the LSND experiment [TBHTB] but 
not confirmed by KARMEN [19] and recently by MiniBooNE [SHIET] . Two well separated 
differences need at least three different neutrino mass eigenstates involved in oscillations. 
Actually the three known neutrino species can be sufficient. At least two u's must be 
massive while, in principle, the third one could still be massless. In the following we will 
assume the simplest picture with three active neutrinos, no sterile neutrinos and CPT 
invariance. The mass eigenstates involved in solar oscillations are mi and m2 and, by 
definition, \m2\ > \mi\, so that Am^^^ = Am^i = 1^2^ — |mip > 0. The atmospheric 
neutrino oscillations involve m^: Am^j^ = |Am|]^| with Amgj^ = jmsp — |mip either 
positive (normal hierarchy) or negative (inverse hierarchy). The present data [221428] are 
compatible with both cases. The degenerate spectrum occurs when the average absolute 
value of the masses is much larger than all mass squared differences: >> |Am^^|. 

With the standard set of notations and definitions [1] the present data are summarised in 
Table 1. 



Quantity 


ref. [251126] 


ref. [271128] 


AmL„ (10-5 eV^) 
Aml,^ (10-3 eV^) 

sin^ 6'i2 

sin^ 6^23 

sin^ 6*13 


7 67+°-i6 

' -0.19 

2 3q+0-ii 

n qi O+0.019 

U.OiZ_Q Q^g 

n 466+°-°'^3 

U.^UU_Q 058 

0.016 ±0.010 


7 65+°-23 

' -"^-0.20 

2 40+0-12 

'^■^^-0.11 
U.OU^_Q 016 

U.UlU_o.011 



Table 1: Fits to neutrino oscillation data. 

Oscillation experiments do not provide information about either the absolute neutrino 
mass scale or the Dirac/Majorana nature of neutrinos. Limits on the mass scale are 
obtained P from the endpoint of the tritium beta decay spectrum, from cosmology (see, 
for example [29]) and from neutrinoless double beta decay {0^/3(3) (for a recent review, see, 
for example [30]). From tritium we have an absolute upper limit of 2.2 eV (at 95% CL.) 
on the mass of electron antineutrino [HI], which, combined with the observed oscillation 
frequencies under the assumption of three CPT-invariant light neutrinos, represents also 
an upper bound on the masses of the other active neutrinos. Complementary information 
on the sum of neutrino masses is also provided by the galaxy power spectrum combined 
with measurements of the cosmic microwave background anisotropics. According to recent 
analyses of the most reliable data [32] \mi\ < 0.60 -^ 0.75 eV (at 95% CL.) depending 
on the retained data (the numbers for the sum have to be divided by 3 in order to obtain a 
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limit on the mass of each neutrino). The discovery of Oz//3/3 decay would be very important 
because it would establish lepton number violation and the Majorana nature of z/'s, and 
provide direct information on the absolute scale of neutrino masses. The present limit from 
Oz//3/3 (with large ambiguities from nuclear matrix elements) is about |mee| < (0.3 -7- 0.8) 
eV [SniEg (see eq. (H). 

After KamLAND [12] , SNO [9] and the upper limits on the absolute value of neutrino 
masses not too much hierarchy in the spectrum of neutrinos is indicated by experiments: 

r = AmL/Am^,^ ~ 1/30 . (1) 

Precisely r = 0.0321qqq5 at 3cr's [251428] . Thus, for a hierarchical spectrum, m2/m^ ~ 
i/r ~ 0.2, which is comparable to the Cabibbo angle Ac ~ 0.22 or to its leptonic analogue 
\fm^Jm'r ~ 0.24. This suggests that the same hierarchy parameter (raised to powers with 
0(\) exponents) may apply for quark, charged lepton and neutrino mass matrices. This 
in turn indicates that, in the absence of some special dynamical reason, we do not expect 
quantities like 6*13 or the deviation of 6*23 from its maximal value to be too small. Indeed 
it would be very important to know how small the mixing angle ^13 is and how close to 
maximal ^23 is. 

Given that neutrino masses are certainly extremely small, it is really difficult from 
the theory point of view to avoid the conclusion that the lepton number L conservation 
is probably violated and that z/'s are Majorana fermions. In this case the smallness of 
neutrino masses can be naturally explained as inversely proportional to the large scale 
where L conservation is violated. If neutrinos are Majorana particles, their masses arise 
from the generic dimension-five non renormalizable operator of the form [33]: 

O. ^ ™%(^ + . (2) 

with H being the ordinary Higgs doublet, /j the SU(2) lepton doublets, 77 a matrix in 
fiavor space, M a large scale of mass and a charge conjugation matrix C between the 
lepton fields is understood. For r/jj C'(l), neutrino masses generated by O5 are of 
the order ^ jM where v ~ 0(100 GeV) is the vacuum expectation value of the 
ordinary Higgs. A particular realization of this effective mass operator is given by the 
see-saw mechanism [331 - 138] , where M derives from the exchange of heavy neutral objects 
of weak isospin or 1. In the simplest case the exchanged particle is the right-handed 
(RH) neutrino z/^ (a gauge singlet fermion here described through its charge conjugate 
field), and the resulting neutrino mass matrix reads (type I see-saw ) [1]: 



m^M-^niD , (3) 



where m£, and M denote the Dirac neutrino mass matrix (defined as z/'^-^m^jz/) and the 
Majorana mass matrix of u'^ (defined as u^'^Mu'^), respectively. As one sees, the light 
neutrino masses are quadratic in the Dirac masses and inversely proportional to the large 
Majorana mass. For mj, ^ y^Amf^ ^ 0.05 eV and m,^ ^ rri^^/M with mo ^ v ^ 
200 GeV we find M ^ 10^^ GeV which indeed is an impressive indication that the scale 
for lepton number violation is close to the grand unified scale Mqut- Thus probably 
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Figure 1: A plot |39| of in cV, the quantity measured in neutrino- less double beta decay, given in 
eq. (U), versus the lightest neutrino mass mi, also in eV. The upper (lower) band is for inverse (normal) 
hierarchy. 

neutrino masses are a probe into the physics near Mqut- This argument, in our opinion, 
strongly discourages models where neutrino masses are generated near the weak scale and 
are suppressed by some special mechanism. 

Oscillation experiments cannot distinguish between Dirac and Majorana neutrinos. 
The detection of neutrino-less double beta decay would provide direct evidence of L non 
conservation, and the Majorana nature of neutrinos. It would also offer a way to possibly 
disentangle the 3 cases of degenerate, normal or inverse hierachy neutrino spectrum. The 
quantity which is bound by experiments on Oi//3/3 is the 11 entry of the v mass matrix, 
which in general, from ttli, = U*mdiagU\ is given by : 

\mee\ = 1(1 - s^g) (miCi2 + maSis) + mge^^'^Sigl , (4) 

where U = Upmns is the mixing matrix, mi^2 are complex masses (including Majorana 
phases) while ms can be taken as real and positive and is the Upmns phase measur- 
able from CP violation in oscillation experiments. Starting from this general formula it 
is simple to derive the bounds for degenerate, inverse hierarchy or normal hierarchy mass 
patterns shown in Fig. 1 [39]. In the next few years a new generation of experiments will 
reach a larger sensitivity on Oi^/3/3 by about an order of magnitude. If these experiments 
will observe a signal this will be compatible with both type of neutrino mass ordering, 
if not, then the normal hierarchy case remains a possibility. Establishing that L is vio- 
lated in particle interactions would also strongly support the possibility that the observed 
baryon asymmetry is generated via leptogenesis, through the out-of-equilibrium, CP and 
L violating decays of the heavy RH neutrinos (see Sect. 10). 

Neutrino mixing is important because it could in principle provide new clues for the 
understanding of the flavor problem. Even more so since neutrino mixing angles show 
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a pattern that is completely different than that of quark mixing: for quarks all mixing 
angles are small, for neutrinos two angles are large (one is even compatible with the 
maximal value) and only the third one is small. For building up theoretical models of 
neutrino mixing one must guess which features of the data are really relevant in order 
to identify the basic principles for the formulation of the model. In particular, it is an 
experimental fact [221428] that within measurement errors the observed neutrino mixing 
matrix [1] is compatible with the so called Tri-Bimaximal (TB) form in eq. ( IT2|l |10] - H3] . 
The best measured neutrino mixing angle 612 is just about la below the TB value tan^ 612 = 
1/2, while the other two angles are well inside the la interval (see table [1]). Thus, one 
possibility is that one takes this coincidence seriously and only considers models where 
TB mixing is automatically a good first approximation. Alternatively one can assume 
that the agreement of the data with TB mixing is accidental. Indeed there are many 
models that fit the data and yet TB mixing does not play any role in their architecture. 
For example, in ref. there is a list of Grand Unified SO(IO) models with parameters 
that can be fitted to the neutrino mixing angles leading to a good agreement with the 
data although most of these models have no built-in relation with TB mixing (see also 
|l5]). Another class of examples is found in ref. [16]. Clearly, for this type of models, 
in most cases different mixing angles could also be accommodated by simply varying the 
fitted values of the parameters. If instead we assume that TB mixing has a real physical 
meaning, then we are led to consider models that naturally produce TB mixing in first 
approximation and only a very special dynamics can lead to this peculiar mixing matrix. 
Discrete non abelian groups (for an introduction see, for example, [ITHSj) naturally emerge 
as suitable flavor symmetries. In fact the TB mixing matrix immediately suggests rotations 
by fixed, discrete angles. It has been found that a broken flavor symmetry based on the 
discrete group A4, (the group of even permutations of 4 elements, which can be seen as 
the invariance group of a rigid regular tetrahedron) appears to be particularly suitable 
to reproduce this specific mixing pattern in leading order (LO). A non exhaustive list of 
papers that discuss the application of to neutrino mixing is given by [l9] - [9^ . The 
choice of this particular discrete group is not unique and, for example, other solutions 
based on alternative discrete flavor symmetries (for example, T' [17 ][55HTU3] . 5*4 [lU4H118j . 
A(27) [TT5l[TT9tiT23] and other groups [T23HT28] ) or continuous flavor symmetries p^9HT3i] 
have also been considered (for other approaches to TB mixing see |135H141] ). but the 
models have a particularly economical and attractive structure, e.g. in terms of group 
representations and of field content. In most of the models ^4 is accompanied by additional 
flavor symmetries, either discrete like Z^v or continuous like U(l), which are necessary to 
eliminate unwanted couplings, to ensure the needed vacuum alignment and to reproduce 
the observed mass hierarchies. Given the set of flavor symmetries and having specified 
the field content, the non leading corrections to TB mixing arising from higher order 
effects can be evaluated in a well defined expansion. In the absence of specific dynamical 
tricks, in a generic model, all the three mixing angles receive corrections of the same 
order of magnitude. Since the experimentally allowed departures of 9i2 from the TB value 
sin^ 6'i2 = 1/3 are small, at most of C(A^) , with the Cabibbo angle, it follows that 
both 6*13 and the deviation of 6*23 from the maximal value are typically expected in these 
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Figure 2: The values of sin^ 6*12 for TB o BM mixing are compared with the data 

models to also be at most of C(A^) 0. A value of 9is ~ 0{Xq) is within the sensitivity of 
the experiments which are now in preparation and will take data in the near future. 

Going back to the possibility that the agreement of the data with TB mixing is acciden- 
tal, we observe that the present data do not exclude a value for 613, i.e. 6*13 ~ 0{Xc), larger 
than generally implied by models with approximate TB mixing. In fact, recent analysis 
of the available data lead to sin^ ^13 = 0.016 ± 0.010 at la [251E6], sin^ ^13 = 0.010^|];[^i? at 
la [271I28], sin^^is = 0.014+°;°^? at la [I42| and sin^ ^13 = 0.010^|];J^^^ at la [US], which 
are compatible with both options. If experimentally it is found that ^13 is near its present 
upper bound, this could be interpreted as an indication that the agreement with the TB 
mixing is accidental. In fact a different empirical observation is that 612 + Xc ~ 7r/4, a rela- 
tion known as quark-lepton complementarity |143yi44] . or similarly 612 + \ fm^Jrnr n /A. 
No compelling model leading, without parameter fixing, to the exact complementarity re- 
lation has been produced so far. Probably the exact complementarity relation is to be 
replaced with something like 6*12 + 0{Xc) ~ 7r/4 or 612 + 0{^Jm^/mT-) ~ 7r/4 (which we 
could call "weak" complementarity |145H164] . If we take any of these complementarity 
relations as a serious hint then a scheme would be relevant where Bimaximal (BM) mixing, 
instead of TB mixing, is the correct first approximation, modified by terms of 0{Xc)- A 
comparison of the TB or BM mixing values with the data on sin^ 612 is shown in Fig. ([2]). 

A very special dynamics is also needed for BM mixing and again discrete symmetry 
groups offer possible solutions. For example, a model |164] based on 6*4, the permutation 
group of 4 elements, naturally leads to BM mixing in LO. This model is built in such a way 
that the dominant corrections to the BM mixing only arise from the charged lepton sector 
at Next-to-the-Leading-Order (NLO) and naturally inherit Ac as the relevant expansion 
parameter. As a result the mixing angles deviate from the BM values by terms of 0{Xc) 
(at most), and weak complementarity holds. A crucial feature of this particular model is 
that only 612 and ^13 are corrected by terms of 0{Xc) while ^23 is unchanged at this order 
(which is essential to make the model agree with the present data). 

Other types of LO approximations for the lepton mixing pattern have been suggested. 
For instance a viable first approximation of the solar mixing angle is also 612 = tan^^(l/</?) 
where = (l + y5)/2 is the golden ratio [I65]. This leads to sin^ ^12 = 1/(1 + if"^) ^ 0.276, 
not far from the allowed range. Another possible connection with the golden ratio has 
been proposed in ref. |166j . In this case cos 6^12 = ^p/^, or sin^ 6^12 = 1/4(3 — (p) ^ 0.345. 
There have been attempts to reproduce these values by exploiting flavor symmetries of 

^By ©(Ap) we mean numericaUy of order A^. As Ac ^ 0.22 a hnear term in Ap with a smahish 
coefficient can easily be ©(A^) 
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Group 


d 


Irr RpDr s 


Presentation 


Ref.'s 




6 


I V 2 


^3 = ^2^ (v4E)2 = 1 




Da 


8 


I1 1/1 2 


= 52 = (yl5)2 = 1 




D7 


14 


1 1' 2 2' 2" 


= 52 = (A5)2 = 1 


fiiil 


Aa 


12 


1 1' 1" 3 

-L, -L , -L , U 


A^ = = (AB)'-^ = 1 


V J 


A^ ~ P5L2(5) 


60 


1, 3, 3', 4, 5 


A^ = B^ = {BAf = 1 


H 


V 


24 


1 r 1" 2 2' 2" 3 


A^ = (ABY = i?2 = 1 = R 


[VIJ 


Sa 


24 


1 1' 2 3 3' 


5M ■ = 52 = (^5)3 ^ ]^ 
TB : A^ = B^ = (5^2)2 = 1 




A (27) ~ X 

1 — * V / 


27 


Xi. ...Xo. 0. 




viii] 




168 


1,3,3,6,7,8 


^3 = 52 = = (P-lA-lPA)^ = 1 


[ix] 


T7 ~ ^7 XI ^3 


21 


l,l',r,3,3 


= B^ = 1, = BA^ 


M 



Table 2: Some small discrete groups used for model building, [i] [T68HI89] : [ii] pl0lfT9T]: [iii] [T92l[T93] : [iv] 

mm-, M [na; h nziinMim]; [vn] piHrniuM) : [viii] ^rmm\ - n nMmi; m ps] . 



icosahedral type |124] , for the first possibility, or of dihedral type [167] for the second case. 

Thus discrete flavor symmetries may play an important role in models of neutrino 
mixing. In particular this is the case if some special patterns indicated by the data 
as possible first approximations, like TB or BM mixing or others, are indeed physically 
relevant. A list of the simplest discrete groups that have been considered for neutrino 
mixing, with some of their properties, is shown in Table 2. In the present review we will 
discuss the formalism and the physics of a non exhaustive list of models of neutrino mixing 
based on discrete symmetries. 



2 Special patterns of neutrino mixing 



Given the PNMS mixing matrix U (we refer the reader to ref. [T] for its general defini- 
tion and parametrisation) , the general form of the neutrino mass matrix, in terms of the 
(complex ^ mass eigenvalues mi, m2, m^, in the basis where charged leptons are diagonal, 
is given by: 

m,y = U*diag{mi,m2,m^)U^ . (5) 

We will present here a number of particularly relevant forms of U and m^, that will be 
important in the following. We start by the most general mass matrix that corresponds 
to 6*13 = and 623 maximal, that is to U given by (in a particular phase convention): 



U 



/ 


C12 


S12 


^ 




Sl2 


C12 


1 




V2 


V2 






S12 


C12 


1 


\ 


V2 


V2 


7i / 



(6) 



^We absorb the Majorana phases in the mass eigenvalues rrii, rather than in the mixing matrix U. The 
dependence on these phases drops in neutrino oscillations. 
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with Ci2 = cos and S12 = sin6'i2. By applying eq. we obtain a matrix of the 
form [T95ti2T0] : 

/x y y\ 

m = \ y z w \ , (7) 
\y w z I 

with complex coefficients x, ?/, z and w. This matrix is the most general one that is 
symmetric under 2-3 (or /i — r) exchange or: 



where ^23 is given by: 



A 



23 




(9) 



The solar mixing angle 9i2 is given by 
sin^ 26 



8 


x*y + y*{w + z)\'^ 


8 


x*y + y*{w + 2;) P + ( w + z 


2 _ 


a; 


2')2 



81/^ 



(x — w — 2)2 + 8?/^ 



(10) 



where the second equality applies to real parameters. Since ^13 = there is no CP violation 
in neutrino oscillations, and the only physical phases are the Majorana ones, accounted 
for by the general case of complex parameters. We restrict here our consideration to 
real parameters. There are four of them in eq. (JTj) which correspond to the three mass 
eigenvalues and one remaining mixing angle, 9i2. Models with yU-r symmetry have been 
extensively studied |195ti211j . 

The particularly important case of TB mixing is obtained when sin^ 2^12 = 8/9 or 
a;-|-?/ = w + 2;§. In this case the matrix takes the form: 



X y y 
m„ = \ y X + V y — V 

y y — V X + V 



'ID 



In fact, in this case, U = Utb is given by 



TB 



\ 



2 1 

-3 7! ° 

1 1 1 

1 1 1 

V V6 V3 V2 ^ 



(12) 



and, from eq. (j5]) one obtains: 



mi^l^^ + m2$2$2 + ^^3*3^ 



(13) 



•^The other solution x — y = w + z gives rise to TB mixing in another phase convention and is physically 
equivalent to x + y = w + z. 
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where 



T 



tj^ ^(1.1.1) 



(14) 



are the respective columns of Utb and are the neutrino mass eigenvalues (mi = x — y, 
m2 = X + 2y and 777,3 = x — y + 2v). It is easy to see that the TB mass matrix in eqs. 
fll3|14p is indeed of the form in eq. (11 ip . All patterns for the neutrino spectrum are 
in principle possible. For a hierarchical spectrum m-^ » m2 » mi, ml ~ Am^ 



''atmi 



ml/ml ~ Am^^i/ Aml^^ and mi could be negligible. But also degenerate masses and 
inverse hierarchy can be reproduced: for example, by taking m^ = —m2 = mi we have a 
degenerate model, while for mi = —m2 and m^ = an inverse hierarchy case is realized 
(stability under renormalization group running (for a review see for example |212] ) strongly 
prefers opposite signs for the first and the second eigenvalue which are related to solar 
oscillations and have the smallest mass squared splitting). 

Note that the mass matrix for TB mixing, in the basis where charged leptons are 
diagonal, as given in eq. (11 ip . can be specified as the most general matrix which is 
invariant under /i — r (or 2-3) symmetry (see eqs. ([H]),(in])) and, in addition, under the 
action of a unitary symmetric matrix Stb (actually S^^ = 1 and [S'tb, ^423] = 0): 



m„ 



SrBmySxB 



A23m^A23 



where Stb is given by: 



S 



TB 




2 

2 



As a last example consider the case of BM where, in addition to 6*13 
maximal, one also has sin^ 2^12 = 1. The BM mixing matrix is given by: 



(15) 

(16) 
and 6*23 



BM 





1 


\ 




"71 


1 


1 


1 


2 


2 




1 


1 


1 


V 2 


2 


71 / 



(17) 



In the basis where charged lepton masses are diagonal, from eq. (|5]), we derive the effective 
neutrino mass matrix in the BM case: 



m„ = mi<l>i<I)f + m2$2$^ + "^3*3$^ 



where 



€-^(0,-1,1) 



'1^ 



(19) 



As we see the most general mass matrix leading to BM mixing is of the form: 



m^, 




(20) 



8 



The resulting matrix can be completely characterized by the requirement of being invariant 
under the action of A23 and also of the unitary, real, symmetric matrix Sbm (satisfying 
= 1 and [Sbm, A23] = 0): 

rriu = SbmIti^Sbm , rri^ = ^ssm^Aas , (21) 

with Sbm given by: 



S. 



BM 





1 

f 



1 


1 


\ 








2 


2 




1 


1 




~2 


2 


/ 



(22) 



V V2 

The nil, niass matrices of the previous examples were all derived in the basis where 
charged leptons are diagonal. It is useful to consider the product = mlrrie, where mg 
is the charged lepton mass matrix (defined as ipp^irieipL), because this product transforms 
as m'^ = U^m^Ue, with Ue the unitary matrix that rotates the left-handed (LH) charged 
lepton fields. The most general diagonal is invariant under a diagonal phase matrix 
with 3 different phase factors: 

mime = T^m^meT (23) 

and conversely a matrix m^m^. satisfying the above requirement is diagonal. If T"^ = 1 the 
matrix T generates a cyclic group Z„. In the simplest case n=3 and we get but n > 3 
is equally possible. Examples are: 

Ttb = I a; I . (24) 



where o;^ = 1, so that T^^ = 1, or 



Tbm=\ -z . (25) 





with T|jv^ = 1. 

We are now in a position to explain the role of finite groups and to formulate the 
general strategy to obtain one of the previous special mass matrices, for example that of 
TB mixing. We must find a group Gj which, for simplicity, must be as small as possible 
but large enough to contain the S and T transformations. A limited number of products 
of S and T close a finite group Gf. Hence the group Gf contains the subgroups Gs and Gt 
generated by monomials in S and T, respectively. We assume that the theory is invariant 
under the spontaneously broken symmetry described by Gf. Then we must arrange a 
breaking of Gj such that, in leading order, Gj is broken down to G5 in the neutrino mass 
sector and down to Gt in the charged lepton mass sector. In a good model this step must 
be realized in a natural way as a consequence of the stated basic principles, and not put 
in by hand. The symmetry under in some cases is also part oi Gf (this the case of 
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S4) and then must be preserved in the neutrino sector along with S by the Gf breaking 
or it could arise as a consequence of a special feature of the Gf breaking (for example, in 
^4 it is obtained by allowing only some transformation properties for the flavons with non 
vanishing VEV's). The explicit example of A4 is discussed in the next section. Note that, 
along the same line, a model with fi — t symmetry can be realized in terms of the group 
5*3 generated by products of A23 and T (see, for example, |189] ). 



3 The A4 group 

is the group of the even permutations of 4 objects. It has 4!/2=12 elements. Geo- 
metrically, it can be seen as the invariance group of a tetrahedron (the odd permutations, 
for example the exchange of two vertices, cannot be obtained by moving a rigid solid). Let 
us denote a generic permutation (1,2,3,4) — (ni, n2, ria, 'n.4) simply by (721712^3714). /I4 
can be generated by two basic permutations S and T given by S* = (4321) and T = (2314). 
One checks immediately that: 

= T^ = {STY = 1 (26) 

This is called a "presentation" of the group. The 12 even permutations belong to 4 
equivalence classes {h and k belong to the same class if there is a in the group such that 
ghg~^ = k) and are generated from S and T as follows: 

CI : / = (1234) (27) 

C2 ■T= (2314), ST = (4132), = (3241), STS = (1423) 

C3 : = (3124), ST^ = (4213), T^^ = (2431), TST = (1342) 

GA ■.S = (4321), T^^T = (3412), TST^ = (2143) 



Note that, except for the identity / which always forms an equivalence class in itself, the 
other classes are according to the powers of T (in C4 5* could as well be seen as ST^). 

The characters of a group xf ^^e defined, for each element g, as the trace of the 
matrix that maps the element in a given representation R. From the invariance of traces 
under similarity transformations it follows that equivalent representations have the same 
characters and that characters have the same value for all elements in an equivalence class. 
Characters satisfy xfXg* = N6^^ , where is the number of transformations in the 
group (A^ = 12 in Also, for each element h, the character of /i in a direct product 
of representations is the product of the characters: Xh^^ = XhXh also is equal to the 
sum of the characters in each representation that appears in the decomposition of i? ® 5". 
In a finite group the squared dimensions of the inequivalent irreducible representations 
add up to A^. The character table of A4 is given in Table 3. From this table one derives 
that A4 has four inequivalent representations: three of dimension one, 1, 1' and 1" and 
one of dimension 3. 

It is immediate to see that the one- dimensional unitary representations are obtained 

by: 

1 ^ = 1 T = 1 (28) 
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Class 






A. 


3 

A 


Ci 


1 


1 


1 


3 


Ca 


1 









C3 


1 









C4 


1 


1 


1 


-1 



Table 3: Characters of Aa 



1' s 
I" s 



1 T 
1 T 



gZ27r/3 
gZ47r/3 



a; 



Note that a; = — 1/2 + iA/3/2 is the cubic root of 1 and satisfies = w*, 1 + a; + = 0. 

The three-dimensional unitary representation, in a basis where the element S* = S" is 
diagonal, is built up from: 

1 \ /O 1 0\ 

S' = \ Q -I , T' = 1 . (29) 
,0 -1/ \1 0/ 

The multiplication rules are as follows: the product of two 3 gives 3x3 = l+l'+l"+3+3 
and 1' X 1' = 1", r X 1" = 1, 1" x 1" = V etc. If 3 ~ (ai, a2, as) is a triplet transforming by 
the matrices in eq. ( I29l) we have that under 5": S"(ai, 02, as)* = (fli, —02, — ^s)* (here the 
upper index t indicates transposition) and under T': T'(ai, a2, as)* = (02, as? c^i)*- Then, 



from two such triplets 3^ ~ (ai,a2,as), 
obtained from their product are: 



(&i5&2,&3) the irreducible representations 



1 = aihi + 0262 + 03^3 
1' = aihi + a;^a2&2 + ^0363 
1" = aihi + u;a2&2 + ^^^03^3 

3 ~ (a2&3,«3^1,«1^2) 

3 ~ (0362,0163,0261) 



(30) 

(31) 
(32) 
(33) 
(34) 



In fact, take for example the expression for 1" = Oi^i + ^0262 + i^^os^s- Under S' it is 
invariant and under T' it goes into 0262 + 1^0363 + a;^ai6i = a;^[ai6i +070262 + 1^^0363] which 
is exactly the transformation corresponding to 1". 

In eq. ( I29l) we have the representation 3 in a basis where S is diagonal. We shall see 
that for our purposes it is convenient to go to a basis where instead it is T that is diagonal. 
This is obtained through the unitary transformation: 



T 



S 



VT'V^ 



VS'V^ 



1 











UJ 















-1 


2 


( 


2 


-1 




2 


2 



2 
2 



(35) 



(36) 
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where: 

y = i= f 1 u \ . (37) 

The matrix V is special in that it is a 3x3 unitary matrix with all entries of unit absolute 
value. It is interesting that this matrix was proposed long ago as a possible mixing matrix 
for neutrinos |213[I214] . We shall see in the following that in the T diagonal basis the 
charged lepton mass matrix (to be precise the matrix m\me) is diagonal. Notice that 
the matrices {S,T) of eqs. ( I35ll36l) coincide with the matrices {Stb,Ttb) of the previous 
section. 

In this basis the product rules of two triplets, {ipi,4'2,4'3) and {ipi,(p2,f3) of A4, ac- 
cording to the multiplication rule 3x3 = 1 + I' + l" + 3 + 3 are different than in the 
S diagonal basis (because for Majorana mass matrices the relevant scalar product is [ab) 
and not (a^b) and are given by: 







V'lV^l + ^2^3 + ^3V^2 


~1 , 








^SV^S + i^lf2 + i^2fl 










1p2^2 + '4'3^1 + 1pl^3 






2V'i(^i 


- ^2^3 


- '^3^2 \ 


/ 'ip2V3 


- 1p3^2 


'2i'3V^3 


- tpl^2 




i'1^2 


- 1p2^1 


2V^2</'2 


- tpl^3 


-'03</'l / 


\ 1p3^l 


- ^1^3 



- 3a . (38) 

In the following we will work in the T diagonal basis, unless otherwise stated. In this 
basis the 12 matrices of the 3-dimensional representation of A4 are given by: 

















C2 : 




























TS 
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C4 : ^=11 2 -1 2 \ ,T^ST 





TST^ 



We can now see why works for TB mixing. In section 2 we have aheady mentioned 
that the most general mass matrix for TB mixing in eq. ( fTTl) . in the basis where charged 
leptons are diagonal, can be specified as one which is invariant under the 2-3 (or /i — r) 
symmetry and under the S unitary transformation, as stated in eq. (fT5|) (note that Stb 
in eqs. (fT5| [T6|) coincides with S in eq. (!36|) ). This observation plays a key role in leading 
to ^4 as a candidate group for TB mixing, because 5* is a matrix of ^4. Instead the 
matrix is not an element of (because the 2-3 exchange is an odd permutation). We 
shall see that in ^4 models the 2-3 symmetry is maintained by imposing that there are no 
flavons transforming as 1' or 1" that break ^4 with two different VEV's (in particular one 
can assume that there are no flavons in the model transforming as 1' or 1"). It is also clear 
that a generic diagonal charged lepton matrix m^me is characterized by the invariance 
under T, or T'^m^rrif.T = m^rrif.. 

The group ^4 has two obvious subgroups: Gg, which is a reflection subgroup generated 
by S and Gt-, which is the group generated by T, which is isomorphic to Z3. If the flavor 
symmetry associated to A^ is broken by the VEV of a triplet (p = {ipi,(p2,f3) of scalar 
fields, there are two interesting breaking pattern. The VEV 

{if) = {vs,vs,vs) (39) 

breaks A4 down to Gs, while 

{ip) = {vT,0,0) (40) 

breaks A4 down to Gt- As we will see, Gs and Gt are the relevant low-energy symmetries 
of the neutrino and the charged-lepton sectors, respectively. Indeed we have already seen 
that the TB mass matrix is invariant under Gs and a diagonal charged lepton mass mime 
is invariant under Gt- 



4 Applying A4 to lepton masses and mixings 

In the lepton sector a typical A4 model works as follows |65]. One assigns leptons to 
the four inequivalent representations of A^. LH lepton doublets I transform as a triplet 
3, while the RH charged leptons e'^, fi'^ and t'^ transform as 1, 1" and 1', respectively. 
Here we consider a see-saw realization, so we also introduce conjugate neutrino fields u'^ 
transforming as a triplet of A4. We adopt a supersymmetric (SUSY) context also to 
make contact with Grand Unification (flavor symmetries are supposed to act near the 
GUT scale). In fact, as well known, SUSY is important in GUT's for offering a solution 
to the hierarchy problem, for improving coupling unification and for making the theory 
compatible with bounds on proton decay. But in models of lepton mixing SUSY also 



13 



helps for obtaining the vacuum ahgnment, because the SUSY constraints are very strong 
and hmit the form of the superpotential very much. Thus SUSY is not necessary but 
it is a plausible and useful ingredient. The flavor symmetry is broken by two triplets 
(fs and ifT and by one or more singlets ^. All these fields are invariant under the SM 
gauge symmetry. Two Higgs doublets /i^^^, invariant under ^4, are also introduced. One 
can obtain the observed hierarchy among me, and m,- by introducing an additional 
U(1)_P7V flavor symmetry [215] under which only the RH lepton sector is charged (recently 
some models were proposed with a different VEV alignment such that the charged lepton 
hierarchies are obtained without introducing a t/(l) symmetry [8TI - I83] ). We recall that 
U{1)fn is a simplest flavor symmetry where particles in different generations are assigned 
(in general) different values of an Abelian charge. Also Higgs flelds may get a non zero 
charge. When the symmetry is spontaneously broken the entries of mass matrices are 
suppressed if there is a charge mismatch and more so if the corresponding mismatch is 
larger. We assign FN-charges 0, q and 2q to r^, ^'^ and e"^, respectively. There is some 
freedom in the choice of q. Here we take q = 2. By assuming that a flavon 6, carrying a 
negative unit of FN charge, acquires a VEV (^)/A = A < 1, the Yukawa couplings become 
fleld dependent quantities ye,^i,T = ye,^l,T{^) and we have 

yr^O{l) , y,^0{\^) , ye^O{\') . (41) 

Had we chosen g = 1, we would have needed {6)/A of order A^, to reproduce the above 
result. The superpotential term for lepton masses, wi is given by: 

wi = 2/ee^(</?T0+y/./^^(v'T/)'+yrr'(v?r/)''+y(i^''0 + (a;A^+^AO(i^V^)+XB((/?sz/V^) + /i.c. + ... 

(42) 

with dots denoting higher dimensional operators that lead to corrections to the LO approx- 
imation. In our notation, the product of 2 triplets (33) transforms as 1, (33)' transforms 
as 1' and (33)" transforms as 1". To keep our formulae compact, we omit to write the 
Higgs and flavon flelds hu^d, and the cut-off scale A. For instance yee'^{ipTl) stands for 
yee^{ipTl)hdO^ / A.^ . The parameters of the superpotential wi are complex, in particular 
those responsible for the heavy neutrino Majorana masses, xa,b- Some terms allowed by 
the symmetry, such as the terms obtained by the exchange ^ (or the term 
iy'^v'^)) are missing in wi. Their absence is crucial and, in each version of models, is mo- 
tivated by additional symmetries (In ref. [56] a natural solution of this problem based on a 
formulation with extra dimensions was discussed; for a similar approach see also [78|[93]). 
In the present version the additional symmetry is Z3. A U{\)r symmetry related to R- 
parity and the presence of driving flelds in the flavon superpotential are common features 
of supersymmetric formulations. Eventually, after the inclusion of N = 1 SUSY breaking 
effects, the U{1)r symmetry will be broken at the low energy scale msusY down to the 
discrete R-parity. Supersymmetry also helps producing and maintaining the hierarchy 
{hu,d) = Vu,d ^ A where A is the cut-off scale of the theory. The flelds in the model and 
their classiflcation under the symmetry are summarized in Table ID 

In this set up it can be shown that the flelds (pT, fs and ^ develop a VEV along the 
directions: 
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I 




















T 


' U 


^0 




3 


1 


1" 


1' 


3 


1 


1 


3 


3 


1 


3 


3 


1 


Zz 




9 




9 




1 


1 


1 






1 






U{1)fn 





4 


2 











-1 




















U{1)r 


1 


1 


1 


1 


1 

















2 


2 


2 



Table 4: Transformation properties of all the fields. 



i'^s) = {vs,vs,vs) 

(0 



u 



(43) 



A crucial part of all serious models is the dynamical generation of this alignment in a 
natural way. We refer to ref. [65] for a proof that the above alignment naturally follows 
from the most general LO superpotential implied by the symmetries of the model. As 
already mentioned, the group A^ has two obvious subgroups: G5, which is a reflection 
subgroup generated by S and Gt, which is the group generated by T, isomorphic to Z^. 
In the basis where S and T are given by eq. f|T6l) . the VEV {^t) = {vt,0,0) breaks A4 
down to Gt, while (ips) = {vs, vs, vs) breaks A4 down to Gs- 

If the alignment in eq. P3|) is realized, at the leading order of the 1/A expansion, the 
mass matrices mi and m,^ for charged leptons and neutrinos correspond to TB mixing. 
The charged lepton mass matrix is diagonal: 



Vt 



The charged fermion masses are given by: 



Ve 

























(44) 



me = VeVd- 



A 



m. 



ViiVd- 



Vt 
A 



rrir 



VrVd 



Vt 
A 



(45) 



where the suppression coming from the breaking of U{1)fn is understood. For example t/e 
stands for UeO^/ A^. In the neutrino sector, after electroweak and ^4 symmetry breaking 
we have Dirac and Majorana masses: 



D 




yvu 



A + 2B/3 -B/3 -B/3 
M= ( -B/3 2B/3 A- B/3 
-B/3 A- B/3 2B/3 



u 



where 



A = 2x, 



B = 2x 



B- 



U 



The eigenvalues of M are 

Mi = {A + B)u 



Mo = Au 



-A + B)u. 



(46) 



(47) 



(45 
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The mass matrix for light neutrinos is mj, = {m^)'^ M ^m^ with eigenvalues 

y'^^l y^^l y'^vl 

The mixing matrix is Utbi eq. f|T2|) . Both normal and inverted hierarchies in the neutrino 
mass spectrum can be realized. It is interesting that A/^ models with the see-saw mechanism 
typically lead to a light neutrino spectrum which satisfies the sum rule (among complex 
masses) : 

112 
— = . 50 

7713 mi m2 

The phases of the complex parameters A and B do not produce any CP violation in 
neutrino oscillations, since ^13 = 0, but are quite important to make the above sum rule 
compatible with the present data on neutrino masses. A detailed discussion of a spectrum 
of this type can be found in refs. [65|[83]. 

Both types of ordering, normal and inverted are allowed and the above sum rule gives 
rise to bounds on the lightest neutrino mass. For normal ordering we have 



mi > ^/^!^ I i_i^r + ... 1 ^ 0.004 eV 



mi < f l + l^r + ... I f^^ 0.006 eV (51) 
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and for the inverted ordering: 



m3 > y (1 - V + ^ 0.017 eV, (52) 

where the dots represent terms with higher powers of r. Notice that for normal ordering 
the neutrino mass spectrum is essentially determined: mi 0.005 eV, mi 0.01 eV and 
m3 ^ 0.05 eV. Also the possible values of |mee| are restricted. For normal hierarchy we 
have ^ 

\mee\ ^ ^^Am^„„ ^ 0.007 eV . (53) 



3^3 



while for inverted hierarchy 



Imeel > w'^I^^ ^ 0.017 eV . (54) 



In a completely general framework, without the restrictions imposed by the fiavor sym- 
metry, |mee| could vanish in the case of normal hierarchy. In this model |mee| is always 
different from zero, though its value for normal hierarchy is probably too small to be 
detected in the next generation of Oi^/3/3 experiments. 

Note that in the charged lepton sector the flavor symmetry A^ is broken by (</9r) down 
to Gt- Actually the above mass terms for charged leptons are the most general allowed 
by the symmetry Gt- At leading order in 1/A, charged lepton masses are diagonal simply 
because there is a low-energy Gt symmetry. In the neutrino sector ^4 is broken down 
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to Gs, though neutrino masses in this model are not the most general ones allowed by 
Gs- The additional property which is needed, the invariance under A^^-, is obtained by 
stipulating that there are no breaking flavons transforming like V and 1". In fact, from 
eq. (l38l) . we see that the expressions for (33)' and (33)" are not 2-3 symmetric. 

At the next level of approximation each term of the superpotential is corrected by 
operators of higher dimension whose contributions are suppressed by at least one power 
of VEV's/A. The corrections to the relevant part of the superpotential determine small 
deviations from the LO VEV alignment configuration. The NLO corrections to mass 
and mixing matrices are obtained by inserting the corrected VEV alignment in the LO 
operators plus the contribution of the new operators evaluated with the unperturbed 
VEV's. The final result is [SSj that, when the NLO corrections are included, TB mixing 
is violated by small terms of the same order for all mixing angles: 

\ + 0{e) (55) 
0{e) 

where e is of order of the typical VEV in units of A. The fact that TB mixing is well 
satisfied by the data sets the restriction e < 0{\y). From the requirement that the 
Yukawa coupling Uj. remains in the perturbative regime, we also get a lower bound on e of 
about 0.01, the exact value depending on tan/3 = Vu/vd and on the largest allowed \yr\- 
Thus we approximately have 

0.01 <e< 0.05 . (56) 

From the see-saw relations in eq. (jlH]), assuming a coupling y of order one, we see that the 
heavy RH neutrino masses are all of order 10^^ GeV, close to the GUT scale. The cut-off 
of the theory can be estimated form eq. (1561) to be close to 10^^ GeV. 

The above results in eqs. fl44ti54p on the lepton mass matrices and the neutrino spec- 
trum refer to the LO approximation. Relations among neutrino masses can be affected 
by NLO corrections but, for e varying in the range of eq. ( l56l) . the bounds ( l5T|52l) do 
not appreciably change (see for example ref. [216j for a numerical study of the deviations 
induced by vacuum misalignment). Also corrections induced by the renormalization group 
evolution of the parameters can modify the above predictions, but only in the case of suf- 
ficiently degenerate mass levels mi and m2 with equal phases, which occurs for inverted 
mass ordering and far from the lower bound fl52p |217] . The expansion parameter e di- 
rectly controls also other observables, such as the CP asymmetries of leptogenesis and the 
rates of lepton flavor violating transitions. This provides an interesting link between the 
physics in the early universe relevant for leptogenesis and the low energy physics accessible 
in current experiments. We will discuss the interplay bewteen discrete flavor symmetries 
and leptogenesis in Sect. 10. 



sin^ 6i2 = 

sin^ 6^23 = 
sin 6'i3 = 
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5 Possible origin of A4 



There is an interesting relation [65] between the model considered so far and the 
modular group. This relation could possibly be relevant to understand the origin of the 
A4 symmetry from a more fundamental layer of the theory. The modular group F is the 
group of linear fractional transformations acting on a complex variable z: 

z — )■ - — , ad — be = 1 , (57) 
cz + d ^ ' 

where a, b, c, d are integers. There are infinite elements in F, but all of them can be 
generated by the two transformations: 

1 . ^ 

s: z ^ , t : z ^ z+1 , (58) 

The transformations s and t in ( l58l) satisfy the relations 

= (stf = 1 (59) 

and, conversely, these relations provide an abstract characterization of the modular group. 
Since the relations ( 1261) are a particular case of the more general constraint ( l59l) . it is clear 
that A4 is a very small subgroup of the modular group and that the representations 
discussed above are also representations of the modular group. In string theory the trans- 
formations (|58|) operate in many different contexts. For instance the role of the complex 
variable z can be played by a field, whose VEV can be related to a physical quantity like 
a compactification radius or a coupling constant. In that case s in eq. f lSS]) represents 
a duality transformation and t in eq. (158|) represent the transformation associated to an 
"axionic" symmetry. 

A different way to understand the dynamical origin of ^4 was recently presented in 
ref. [nn] where it is shown that the /I4 symmetry can be simply obtained by orbifolding 
starting from a model in 6 dimensions (6D). In this approach A4 appears as the remnant 
of the reduction from 6D to 4D space-time symmetry induced by the special orbifolding 
adopted. This approach suggests a deep relation between flavor symmetry in 4D and 
space-time symmetry in extra dimensions. 

The orbifolding is defined as follows. We consider a quantum field theory in 6 di- 
mensions, with two extra dimensions compactified on an orbifold T"^ jZ^- We denote by 
2; = X5 + ixg the complex coordinate describing the extra space. The torus is defined 
by identifying in the complex plane the points related by 



z -f z^\ 



1 (60) 



I — 

z ^ z + ^ 7 = e3 , 

where our length unit, 2tiR, has been set to 1 for the time being. The parity Z2 is defined 
by 

-z (61) 

and the orbifold T"^ /Z2 can be represented by the fundamental region given by the trian- 
gle with vertices 0,1,7, see Fig. 3. The orbifold has four fixed points, (zi, 2:2, 2:3, 24) = 
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5Cg 



7 



7 + 1 




Figure 3: Orbifold T2/Z2. The regions with the same numbers are identified with each other. The four 
triangles bounded by solid lines form the fundamental region, where also the edges with the same letters 
are identified. The orbifold T2/Z2 is exactly a regular tetrahedron with 6 edges a,b,c,d,e, f and four 
vertices zi, Z2, 23, Z4, corresponding to the four fixed points of the orbifold. 

(1/2, (1 + 7)/2, 7/2, 0). The fixed point 24 is also represented by the vertices 1 and 7. In 
the orbifold, the segments labelled by a in Fig. 1, (0, 1/2) and (1, 1/2), are identified and 
similarly for those labelled by b, (1, (1 + 7)/2) and (7, (1 + 7)/2), and those labelled by c, 
(0,7/2), (7,7/2). Therefore the orbifold is a regular tetrahedron with vertices at the four 
fixed points. 

The symmetry of the uncompactified 6D space time is broken by compactification. 
Here we assume that, before compactification, the space-time symmetry coincides with the 
product of 6D translations and 6D proper Lorentz transformations. The compactification 
breaks part of this symmetry. However, due to the special geometry of our orbifold, a 
discrete subgroup of rotations and translations in the extra space is left unbroken. This 
group can be generated by two transformations: 



thus generating the group A4. From the previous equations we immediately verify that 
S and T satisfy the characteristic relations obeyed by the generators of A^: = = 
{STY = 1- These relations are actually satisfied not only at the fixed points, but on the 
whole orbifold, as can be easily checked from the general definitions of S and T in eq. 
(!62|) . with the help of the orbifold defining rules in eqs. (160|) and (!6T]) . 

We can exploit this particular geometry of the internal space to build a model with 
^4 flavor symmetry. There are 4D branes at the four fixed points of the orbifolding and 



S : z ^ z + ^ 
T : z ^ uz 




(62) 



Indeed S and T induce even permutations of the four fixed points: 



S : {zi, Z2, Z3, Zi) — )■ (2:4, 2:3, Z2-, Zi) 

T: (^1, ^2; ^35 2:4) — )■ (2:2; 2:3, ^1, 2:4) 



(63) 
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the tetrahedral symmetry of A4 connects these branes. The standard model fields have 
components on the fixed point branes while the scalar fields necessary for the ^4 breaking 
are in the bulk. Each brane field, either a triplet or a singlet, has components on all of the 
four fixed points (in particular all components are equal for a singlet) but the interactions 
are local, i.e. all vertices involve products of field components at the same space-time 
point. In the low-energy limit this model coincides with the one illustrated in the previous 
section. Unfortunately in such a limit the 6D construction does not provide additional 
constraints or predictions. 

This construction can be embedded in a SU(5) GUT |218j . Other discrete groups can 
arise from the compactification of two extra dimensions on orbifolds and the possibilities 
have been classified in |219y22Uj within a field theory approach. In string theory the fiavor 
symmetry can be larger than the isometry of the compact space. For instance in heterotic 
orbifold models the orbifold geometry combines with the space group selection rules of 
the string, as shown in |221] . Discrete fiavor symmetries from magnetized/intersecting 
D-branes are discussed in |222j . Discrete symmetries can also arise from the spontaneous 
breaking of continuous ones. Such a possibility has been discussed in ref. |133[I134] . 

6 Alternative routes to TB mixing 

While ^4 is the minimal fiavor group leading to TB mixing, alternative fiavor groups 
have been studied in the literature and can lead to interesting variants with some specific 
features. 

Recently, in ref. |223] . the claim was made that, in order to obtain the TB mixing 
"without fine tuning", the finite group must be 5*4 or a larger group containing 5*4. For 
us this claim is not well grounded being based on an abstract mathematical criterium 
for a natural model (see also |224] ). For us a physical field theory model is natural if 
the interesting results are obtained from the most general lagrangian compatible with the 
stated symmetry and the specified representation content for the fiavons. For example, 
we obtain from A4 (which is a subgroup of S4) a natural (in our sense) model for the 
TB mixing by simply not including symmetry breaking fiavons transforming like the 1' 
and the 1" representations of A4. This limitation on the transformation properties of the 
fiavons is not allowed by the rules specified in ref. |223] which demand that the symmetry 
breaking is induced by all possible kinds of fiavons (note that, according to this criterium, 
the SM of electroweak interactions would not be natural because only Higgs doublets are 
introduced!). Rather, for naturalness we also require that additional physical properties 
like the VEV alignment or the hierarchy of charged lepton masses also follow from the 
assumed symmetry and are not obtained by fine tuning parameters: for this actually /I4 
can be more effective than 5*4 because it possesses three different singlet representations 
1, 1' and 1". 

Models of neutrino mixing based on 5*4 have in fact been studied |104H118] . The group 
of the permutations of 4 objects S4 has 24 elements and 5 equivalence classes (the character 
table is given in Table E]) that correspond to 5 inequivalent irreducible representations, two 
singlets, one doublet, two triplets: li, I2, 2, 3i and 32 (see Table 2). Note that the squares 
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of the dimensions af all these representations add up to 24. 



Class 




X(l2) 


X(2) 


X(3i) 


X(32) 


Ci 


1 


1 


2 


3 


3 


C2 


1 


1 


2 


-1 


-1 


C3 


1 


-1 





1 


-1 




1 


1 


-1 








C5 


1 


-1 





-1 


1 



Table 5: Characters of 6*4 



For models of TB mixing, one starts from the 5*4 presentation = = [BA^Y = 1 
and identifies, up to a similarity transformation, B^ = S and A = T, where S and T are 
given in eqs. ( IT6l IMl) . In this presentation one obtains a realisation of the 3-dimensional 
representation of 5*4 where the 5* and A23 matrices in eq. f|T5|) that leave invariant the 
TB form of my in eq. f lTTj) as well as the matrix T in eq. of invariance for mime, 

all explicitly appear |110] . In 6*4 the V and 1" of ^4 are collected in a doublet. When 
the VEV of the doublet flavon is aligned along the Gs preserving direction the resulting 
couplings are 2-3 symmetric as needed. In ^4 the 2-3 symmetry is only achieved if the 1' 
and 1" VEV's are identical (which is the 5*4 prediction). As discussed in ref. |110] . in the 
leptonic sector the main difference between A4 and 5*4 is that, while in the typical versions 
of ^4 the most general neutrino mass matrix depends on 2 complex parameters (related to 
the couplings of the singlet and triplet flavons) , in 6*4 it depends on 3 complex parameters 
(because the doublet is present in addition to singlet and triplet flavons). 

Other flavor groups have been considered for models of TB mixing. Some of them in- 
clude S'4 as a subgroup, like PSL-^iy) (the smallest group with complex triplet representa- 
tions) [T25HT27] . while others, like A(27) (which is a discrete subgroup of SU{?,)) [Tl9tiT23] 
or Zj XI Z3 |128] . have no direct relation to S'4 |225] . In Sect. 8 we will consider S'4 again 
in the different context of BM with large corrections from the lepton sector. 

A different approach to TB mixing has been proposed and developed in different ver- 
sions by S. King and collaborators over the last few years |129H132]I225] . The starting point 
is the decomposition of the neutrino mass matrix given in eqs. ( 1T31IT41) corresponding to 
exact TB mixing in the diagonal charged lepton basis: 

mj, = mi$i$f + m2$2'^>2' + ■'^s'^s'^s (64) 

where $f = ^(2, -1, -1), $1" = -^(1, 1, 1), $1" = ^(0, -1,1), are the respective columns 
of Utb and mi are the neutrino mass eigenvalues. Such decomposition is purely kinematical 
and does not possess any dynamical or symmetry content. In the King models the idea 
is that the three columns of Utb are promoted to fiavon fields whose VEVs break the 
family symmetry, with the particular vacuum alignments along the directions Eq. 

directly arises in the see-saw mechanism, = m'^M'^mo, written in the diagonal 
RH neutrino mass basis, M = diag(Mi, M2, M3) when the Dirac mass matrix is given 
by mjj = {vi^i,V2^2iV3^3), where Vi are mass parameters describing the size of the 
VEVs. In this way, to each RH neutrino eigenvalue Mj, a particular light neutrino mass 
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rrii is associated. In the case of a strong neutrino hierarchy this idea can be combined 
with the framework of "Sequential Dominance", where the hghtest RH neutrino, with 
its symmetry properties fixes the heaviest hght neutrino and so on. For no pronounced 
hierarchy the correspondence between Mi and rrii can still hold and one talks of "Form 
Dominance" [226j . In these models the underlying family symmetry of the Lagrangian Gj 
is completely broken by the the combined action of the $j VEV's, and the flavor symmetry 
of the neutrino mass matrix emerges entirely as an accidental residual symmetry of the 
quadratic form of eq. (16^ [225] . The symmetry Gf plays a less direct role and the name 
"Indirect Models" is used by the authors. 

7 Extension to quarks and GUT's 

Much attention has been devoted to the question whether models with TB mixing in 
the neutrino sector can be suitably extended to also successfully describe the observed 
pattern of quark mixings and masses and whether this more complete framework can be 
made compatible with (supersymmetric) SU(5) or SO(IO) Grand Unification. For models 
with approximate TB mixing in the leptonic sector we first consider the extension to 
quarks without Grand Unification and then the more ambitious task of building grand 
unified models. In GUT models based on SU{5) ® Gf or SO{10) ® Gf 0, where Gf is a 
flavor group, clearly all fields in a whole representation of SU{5) or 50(10) must have 
the same transformation properties under Gf. This poses a strong constraint on the way 
quarks and leptons have to transform under Gf. 

7.1 Extension to quarks without GUT's 

The simplest attempts of directly extending models based on to quarks have not 
been satisfactory. At first sight the most appeahng possibility is to adopt for quarks 
the same classification scheme under that one has used for leptons (see, for example, 
|l9] - [5T| |65] ) . Thus one tentatively assumes that LH quark doublets Q transform as a 
triplet 3, while the RH quarks {u'^,d'^), {c'^,s'^) and {t'^,b'^) transform as 1, 1" and 1', 
respectively. This leads to Vu = Vd and to the identity matrix for Vckm = KJ^rf 
lowest approximation. This at first appears as very promising: a LO approximation where 
neutrino mixing is TB and Vckm = 1 is a very good starting point. But there are some 
problems. First, the corrections to Vckm = 1 turn out to be strongly constrained by 
the leptonic sector, because lepton mixing angles are very close to the TB values, and, in 
the simplest models, this constraint leads to a too small Vus (i-e. the Cabibbo angle is 
rather large in comparison to the allowed shifts from the TB mixing angles) [65]. Also in 
these models, the quark classification which leads to Vckm = 1 is not compatible with A4 
commuting with SU(5). An additional consequence of the above assignment is that the 
top quark mass would arise from a non-renormalizable dimension five operator. In that 
case, to reproduce the top mass, we need to compensate the cutoff suppression by some 

■^Thc Pati-Salam group SU{4) (g) SU{2) (g) SU{2) has also been considered, for example in P?fl[^ 
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extra dynamical mechanism. Alternatively, we have to introduce a separate symmetry 
breaking parameter for the quark sector, sufficiently close to the cutoff scale. 

Due to this, larger discrete groups have been considered for the description of quarks. A 
particularly appealing set of models is based on the discrete group T', the double covering 
group of A4 |171[nSHID3]- As we see in Table 2 the representations of T' are those of 
A4 plus three independent doublets 2, 2' and 2". The doublets are interesting for the 
classification of the ffist two generations of quarks |229H232] . For example, in ref. |102j a 
viable description was obtained, i.e. in the leptonic sector the predictions of the A4 model 
are reproduced, while the T' symmetry plays an essential role for reproducing the pattern 
of quark mixing. But, again, the classification adopted in this model is not compatible 
with Grand Unification. 

7.2 Extension to quarks within GUT's 

As a result, the group A4 was considered by many authors to be too limited to also 
describe quarks and to lead to a grand unified description. It has been recently shown [79] 
that this negative attitude is not justified and that it is actually possible to construct a 
viable model based on ^4 which leads to a grand unified theory (GUT) of quarks and 
leptons with TB mixing for leptons and with quark (and charged lepton) masses and 
mixings compatible with experiment. At the same time this model offers an example of an 
extra dimensional SU(5) GUT in which a description of all fermion masses and mixings 
is accomplished. The formulation of SU(5) in extra dimensions has the usual advantages 
of avoiding large Higgs representations to break SU(5) and of solving the doublet-triplet 
splitting problem. The choice of the transformation properties of the two Higgses 
and f/5 has a special role in this model. They are chosen to transform as two different A^ 
singlets 1 and 1'. As a consequence, mass terms for the Higgs colour triplets are not directly 
allowed and their masses are introduced by orbifolding, a la Kawamura |233H235] . In this 
model, proton decay is dominated by gauge vector boson exchange giving rise to dimension 
six operators, while the usual contribution of dimension five operators is forbidden by the 
selection rules of the model. Given the large Mqut scale of SUSY models and the relatively 
huge theoretical uncertainties, the decay rate is within the present experimental limits. A 
see-saw realization in terms of an /I4 triplet of RH neutrinos v'^ ensures the correct ratio 
of light neutrino masses with respect to the GUT scale. In this model extra dimensional 
effects directly contribute to determine the flavor pattern, in that the two lightest tenplets 
Ti and T2 are in the bulk (with a doubling Tj and T/, z = 1, 2 to ensure the correct zero 
mode spectrum), whereas the pentaplets F and T3 are on the brane. The hierarchy of 
quark and charged lepton masses and of quark mixings is determined by a combination 
of extra dimensional suppression factors and of U{1)fn charges, both of which only apply 
to the ffist two generations, while the neutrino mixing angles derive from ^4 in the usual 
way. If the extra dimensional suppression factors and the U{1)fn charges are switched off, 
only the third generation masses of quarks and charged leptons survive. Thus the charged 
fermion mass matrices are nearly empty in this limit (not much of A^ effects remain) and 
the quark mixing angles are determined by the small corrections induced by those effects. 
The model is natural, since most of the small parameters in the observed pattern of masses 
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and mixings as well as the necessary vacuum alignment are justified by the symmetries of 
the model. However, in this case, like in all models based on f/(l)^jv; the number of 0{1) 
parameters is larger than the number of measurable quantities, so that in the quark sector 
the model can only account for the orders of magnitude (measured in terms of powers of 
an expansion parameter) and not for the exact values of mass ratios and mixing angles. 
A moderate fine tuning is only needed to enhance the Cabibbo mixing angle between the 
first two generations, which would generically be of 0{Xq). 

The problem of constructing GUT models based on SU{5) (8> Gj or 5*0(10) ® Gf with 
approximate TB mixing in the leptonic sector has been considered by many authors (see, 
for example [791IM1 - I9 11194] . based on A4). In our opinion most of the models are incomplete 
(for example, the crucial issue of VEV alignment is not really treated in depth as it should) 
and/or involve a number of unjustified, ad hoc fine tuning of parameters. An interesting 
model based on SU{5) ® T' is discussed in ref. |103] . In this model the SU{5) tenplets 
T3 and Ta (a = 1, 2) of the third and of the first two generations are classified as 1 and 
2 of T', respectively, while the SU{5) pentaplets are in a 3 of T'. This model provides a 
good description of fermion masses and mixings and appears simpler than the model in 
ref. [79], which is also based on SU{5). However, the model of ref. |103] is fine tuned. In 
fact one does not understand how it is possible that, for example, the electron and the 
muon masses can come out so widely different as observed, given that in this model their 
left and right components separately transform in an identical way under T'. The reason 
is that in the second term of eq. 5 of ref. |103] . only one of three possible contractions 
has been taken into account. If the missing ones, which are also allowed by the assumed 
symmetry properties, are included with generic coefficients, one in fact finds that the e 
and fi masses are of the same order in the absence of fine tuning. Given that the expansion 
parameter in the model is of 0{Xc) the fine tuning which is needed is large. One possible 
way out would be to invoke some ultraviolet completion of the model where particular 
heavy field exchanges could justify the presence of only the desired couplings after the 
heavy fields are integrated out. Also, in the model of ref. |103] there is no discussion of the 
origin of the required vacuum alignment. Recently some GUT models based on SU(5)xS'4 
have appeared |117]I118] . Also in these models the first two generation fermions are in the 
same 5*4 representations (either a doublet, for tenplets, or a triplet, for pentaplets). In the 
absence of an additional principle the electron and muon mass should naturally be of the 
same order. In ref. |118] the vanishing of the electron mass at LO is obtained by the ad 
hoc choice of one particular minimum of the scalar potential among a continuous family 
of degenerate solutions (see their eqs. (70-71)). In the case of ref. |117] the problem is 
solved by introducing new heavy particles with suitable interactions that, once integrated 
out, produce the desired structure for the mass matrix. 

As for the models based on S'O(IO) (8) Gf we select two recent examples with Gf = S4 
[Iiailll] and Gf = PSL2{7) [HZ]. Clearly the case of ^O(IO) is even more difficult 
than that of 5'f/(5) because the neutrino sector is tightly related to that of quarks and 
charged leptons as all belong to the 16 of 5*0(10) (for a general analysis of 50(10) ® ^4 
see [91]). The strategy adopted in refs. |113[I114[I127] as well as in other 5*0(10) models, 
is as follows. One considers renormalisable fermion mass terms with Higgs multiplets of 
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the SO{10) 10 {h terms) and 126 (/ terms) representations. The Majorana neutrino mass 
matrix arises from the 126. One assumes that the dominant contribution to the Dirac 
masses of fermions is from the h terms with small corrections from the / terms. In first 
approximation the h contribution is a matrix of rank 1 with only the third generation mass 
being non vanishing. The light fermion masses and the quark mixings then arise from the 
/ terms (and from some possible extra terms). The third family dominance is obtained 
by a term with a double flavon factor in ref. |113[I114] (based on 5*4) which then makes 
particularly difficult to keep the corrective / terms small (for this fine tuning is needed 
or a suitable ultraviolet completion). In ref. |127j the dominant h terms are induced by a 
single PSL2{7) sextet fiavon (the existence of complex 3-dim and of 6-dim representations 
is the peculiarity of PSL2{7)). In both models in the neutrino sector one has a sum of 
type I and II see-saw contributions of the form: 

m,, = fvL - nr^-^niD (65) 

where the first term is from the exchange of a triplet Higgs with VEV proportional to Vl 
while the second term is from type I see-saw with RH mass proportional to vr. One must 
assume that the first term is dominant and the second is negligible. Then the leading 
approximation for the fermion Dirac masses is from the h terms and for neutrino masses 
from the / terms. The / terms are diagonalized by the TB mixing unitary matrix. In this 
way the connection between quarks and neutrinos is relaxed and a completely different 
pattern of mixing can be realized in the two sectors. Clearly for the fvi dominance in eq. 
(l65l) one needs vl » v^/vr with v ^ h/f . This needs widely different scales for vl and 
vr in the model and much of the description of the corresponding dynamics, along the lines 
of refs. |236[I237] . remains to be studied in detail. In both of these models the discussion 
of the alignment is not satisfactory. In particular in ref. |127] it is only proven that the 
arbitrary coefficients appearing in the most general allowed superpotential can be fitted 
to lead to the required ratios of components in the VEV's (while for a natural model one 
would require that the alignment automatically follows in a whole region of the parameter 
space). In conclusion, in our opinion, the problem of constructing a satisfactory natural 
model based on 5*0(10) with built-in TB mixing at the LO approximation, remains open. 

8 The S4 group and BM mixing 

If one takes the alternative view that the agreement with TB mixing is accidental and is 
rather oriented to consider "weak" complementarity as a more attractive guiding principle, 
then a better starting point could be BM mixing. In the BM scheme tan^6'i2 = 1, to be 
compared with the latest experimental determination: tan^ 612 = 0.45±0.04 (at la) |22] - [28] 
(see fig. 2), so that a rather large non leading correction is needed, as already mentioned. A 
discrete group approach can also work for BM mixing. We now summarise a model |164] 
based on S4 that leads to BM mixing in first approximation while the agreement with the 
data is restored by large NLO corrections that arise from the charged lepton sector. 

The group S4 is particularly suitable for reproducing BM mixing in LO because the 
unitary matrices Sbm, given in eq. (l22l) . and Tbm, given in eq. (l25ll . directly provide a 
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presentation of S4. We recall that Sbm leaves invariant the most general mass matrix for 
BM mixing in the charged lepton diagonal basis, eq. fl20|) . while Tbm leaves invariant the 
most general diagonal matrix mime for charged leptons (see eqs. fl2T| [23]) ). In fact, from 
Table 2, we see that a possible presentation of 5*4 is given by: 

A^ = B^ = {ABf = 1 (66) 

In terms of 3x3 matrices, we can make the identifications A = Tbm and B = Sbm and eqs. 
fl66p are satisfied. As was the case for the ^4 models, again in this model the invariance 
under A23, which is also necessary to specify BM mixing according to eq. fl?Il) . arises 
accidentally as a consequence of the specific field content and is limited to the contribution 
of the dominant terms to the neutrino mass matrix. 

In the model the 3 generations of LH lepton doublets / and of RH neutrinos u'^ to two 
triplets 3, while the RH charged leptons e*^, fi'^ and t'^ transform as 1, 1' and 1, respectively. 
The 5*4 symmetry is then broken by suitable triplet flavons. All the fiavon fields are singlets 
under the Standard Model gauge group. Additional sjTiimetries are needed, as usual, to 
prevent unwanted couplings and to obtain a natural hierarchy among me, m^ and m^. 
The complete flavor symmetry of the model is 6*4 x Z4 x U{1)fn- A fiavon 6, carrying a 
negative unit of the ?7(l)i?Ar charge F, acquires a vacuum expectation value (VEV) and 
breaks U{1)fn- A supersymmetric context is adopted, so that two Higgs doublets hu,d, 
invariant under 6*4, are present in the model as well as the U{1)r symmetry related to 
R-parity and the driving fields in the fiavon superpotential. Supersymmetry also helps 
producing and maintaining the hierarchy (hu^d) = Vu,d ^ A where A is the cut-off scale of 
the theory. 

The fields in the model and their classification under the symmetry are summarized in 
Table [6l The fields ipf, , ■CjJ and (p'l are the driving fields. The complete superpotential 
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Table 6: Transformation properties of all the fields. 



can be written as w = wi + + Wd- The Wd term is responsible for the alignment. It was 
discussed in ref. |164] and this discussion will not be repeated here. The terms wi and Wu 
determine the lepton mass matrices (we indicate with (. . .) the singlet 1, with (. . .)' the 
singlet 1' and with (. . .)v (V = 2, 3, 3') the representation V) 

(1) n2 (2) n2 (3) n2 
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(67) 
(68) 



where a and h are complex coefficients. Again, to keep our formulae compact, we omit 
to write the Higgs fields hu4- For instance yrT^il^ii) / K stands for yT-T'^{l(f)i)hd/ A, yi^uH) 
stands for y{i'H)hu- The powers of the cutoff A also take into account the presence of 
the omitted Higgs fields. Note that the parameters M, M^, and M'^ defined above 
are dimensionless. In the above expression for the superpotential w, only the lowest order 
operators in an expansion in powers of 1/A are explicitly shown. Dots stand for higher 
dimensional operators that will be discussed later on. The stated symmetries ensure that, 
for the leading terms, the fiavons that appear in wi cannot contribute to and viceversa. 

The potential corresponding to Wd possesses an isolated minimum for the following 
VEV configuration: 



) I 1 1 4 iXl) 
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A I A 



where the factors A, B, C, D should obey to the relations: 

v^/iA2 + + f^AB = (71) 



92 2gigi 
Similarly, the Froggatt-Nielsen fiavon 6 gets a VEV, determined by the D-term associated 
to the local U{1)fn symmetry, and it is denoted by 

f = *. (T3) 

With this VEV's configuration, the charged lepton mass matrix is diagonal 

/ {yi'^B'~yi'^A^ + yi'^AB)t'' \ 
mi= y^Bt 

V yrA J 

so that at LO there is no contribution to the Upmns mixing matrix from the diagonal- 
ization of charged lepton masses. In the neutrino sector for the Dirac and RH Majorana 
matrices we have 

1 \ / 2M + 2aD -2bC -2bC \ 

= { 1 \yvu Mat = -2bC 2M + 2aD A . 

10/ V 2M + 2aD / 

(75) 



Vd (74) 
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The matrix M^r can be diagonalized by the BM mixing matrix Ubm, which represents the 
full lepton mixing at the LO, and the eigenvalues are 

Ml = 2\M + aD - V2bC\A = 2\M + aD + V2bC\A M3 = 2|M + aL>|A . (76) 

After see-saw, since the Dirac neutrino mass matrix commutes with and its square 
is a matrix proportional to unity, the light neutrino Majorana mass matrix, given by the 
see-saw relation mj, = (m^)"^M^^m^, is also diagonalized by the BM mixing matrix and 
the eigenvalues are 

= := r \fn2\ = := r ^".3 = — ; tt • 

2\M + aD - ^/2bC\ A 2\M + aD + V^bC] A 2\M + aD\ A 

(77) 

The light neutrino mass matrix depends on only 2 effective parameters, at LO, indeed the 
terms M and aD enter the mass matrix in the combination F = M + aD. The coefficients 
ye \ Ufj^, Ur, y, a and b are all expected to be of 0(1). A priori M could be of 0{1), 
corresponding to a RH neutrino Majorana mass of 0{A), but, actually, it must be of the 
same order as C and D. In the context of a grand unified theory this would correspond 
to the requirement that M is of 0{Mgut) rather than of 0{Mpianck)- 

We expect a common order of magnitude for the VEV's (scaled by the cutoff A): 

Ar^ B r^v , C ^ D ^v' . (78) 

However, due to the different minimization conditions that determine {A^B) and {C,D), 
we may tolerate a moderate hierarchy between v and v'. Similarly the order of magnitude 
of t is in principle unrelated to those of v and v'. It is possible to estimate the values of v 
and t by looking at the mass ratios of charged leptons (while v' only enters in the neutrino 
sector): and the result is that t ~ 0.06 and v ~ 0.08 (modulo coefficients of 0{1)). 

So far we have shown that, at LO, we have diagonal and hierarchical charged leptons 
together with the exact BM mixing for neutrinos. It is clear that substantial NLO correc- 
tions are needed to bring the model to agree with the data on 612- A crucial feature of the 
model is that the neutrino sector fiavons (pi^ and are invariant under Z4 which is not 
the case for the charged lepton sector fiavons 0/ and xi- The consequence is that 0j, and 
^i, can contribute at NLO to the corrections in the charged lepton sector, while at NLO (pi 
and xi cannot modify the neutrino sector couplings. As a results the dominant corrections 
to the BM mixing matrix only occur at NLO through the diagonalization of the charged 
leptons. In fact, at NLO the neutrino mass matrix is still diagonalized by Ubm but the 
mass matrix of charged leptons is no more diagonal. Including these additional terms from 
the diagonalization of charged leptons the Upmns niatrix can be written as 

UpMNS = UjUsM , (79) 

and therefore the corrections from Ui affect the neutrino mixing angles at NLO according 
to 

sin^ 023 = ^ (80) 
smeis = ^{Vi2-Vi3)v' . 



28 



where the coefficients Vij arise from Ui. By comparing these expressions with the current 
experimental values of the mixing angles in Table 1, we see that, to correctly reproduce 
6*12 we need a parameter v' of the order of the Cabibbo angle Ac*. Moreover, barring 
cancellations of/among some the Vij coefficients, also ^13 is corrected by a similar amount, 
while 623 is unaffected at the NLO. A salient feature of this model is that, at NLO accuracy, 
the large corrections of 0{Xc) only apply to Ou and ^13 while 6*23 is unchanged at this 
order. As a correction of 0{Xc) to 623 is hardly compatible with the present data (see 
Table 1) this feature is very crucial for the phenomenological success of this model. It is 
easy to see that this essential property depends on the selection in the neutrino sector of 
fiavons C,u and (p^, that transform as 1 and 3 of 6*4, respectively. If, for example, the singlet 
C,u is replaced by a doublet ipi, (and correspondingly the singlet driving field is replaced 
by a doublet all other quantum numbers being the same, one can construct a variant 
of the model along similar lines, but, in this case, all the 3 mixing angles are corrected 
by terms of the same order. This confirms that a particular set of 6*4 breaking fiavons is 
needed in order to preserve 6'23 from taking as large corrections as the other two mixing 
angles. 

All this discussion applies at the NLO and we expect that at the NNLO the value of 
$23 will eventually be modified with deviations of about (!?(A^). The next generation of 
experiments, in particular those exploiting a high intensity neutrino beam, will probably 
reduce the experimental error on 623 and the sensitivity on ^13 to few degrees. All quanti- 
tative estimates are clearly affected by large uncertainties due to the presence of unknown 
parameters of order one, but in this model a value of 613 much smaller than the present 
upper bound would be unnatural. If in the forthcoming generation of experiments no 
significant deviations from zero of ^13 will be detected, this construction will be strongly 
disfavoured. 

9 Lepton flavor violation 

Neutrino oscillations provide evidence of fiavor conversion in the lepton sector. This 
indicates that lepton fiavor violation (LFV) might take place, at least at some level, also in 
other processes such as those involving charged leptons. Flavor violating decays of charged 
leptons, strictly forbidden in the SM, are indeed allowed as soon as neutrino mass terms 
are considered. If neutrino masses are the only source of LFV, the effects are too small to 
be detected, but in most extensions of the SM where new particles and new interactions 
with a characteristic scale A^p are included, the presence of new sources of flavor violation, 
in both quark and lepton sectors, is a generic feature. The scale A^p can be much smaller 
than the cut-off scale A introduced before. Indeed there are several indications suggesting 
new physics at the TeV scale, such as a successful gauge coupling uniflcation, viable 
solutions to the hierarchy problem and realistic dark matter candidates. In a low-energy 
description, the associated effects can be parametrized by higher-dimensional operators. 
The dominant terms are represented by dimension six operators, suppressed by two powers 
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of Atvp: 




A2 



2 — e'^iH^ a^'^ Ffj_^Zijlj + — ^ — [4 — f ermion operators] + h.c. (81) 




where e is the electric charge and Zij denotes an adimensional complex matrix with indices 
in flavor space. If the underlying theory is weakly interacting with a typical coupling 
constant gNP and predicts new particles of mass uinp we expect A^p f» AnmNp/ Qnp- 
The present bounds on the branching ratios |238] of the rare charged lepton decays set 
stringent limits on combinations of the scale A^rp and the coefficients of the involved 
operators. For instance, from BR{^ 67) < 1.2 x 10"^^ |239[I240] we get |Z^e| < 10"*^ x 
[A7vp(TeV)/l TeV]^. Typically, for coefficients of order one, the existing bounds require 
a large scale A^vp, several orders of magnitude larger than the TeV scale. Conversely, 
to allow for new physics close to the TeV scale, coefficients much smaller than one are 
required, which may indicate the effect of a flavor symmetry. 

In theories with a flavor symmetry group Gf spontaneously broken by a set of small 
parameters e, the coefficients of the effective lagrangian in eq. (18 ip become functions of 
e. The low-energy Lagrangian of eq. (ISTj) is derived from the theory defined close to the 
cut-off scale A, where all operators are invariant under Gf thanks to their dependence on 
the flavon multiplets. Below the flavor symmetry breaking scale the flavons are replaced 
by their VEVs, which enter the coefficients of L^ff through the dimensionless combination 
e ~ VEV/ A. Exploiting the smallness of the parameters e we can keep in >Ce// the first 
few terms of a power series expansion. For instance: 



Notice that the same symmetry breaking parameters that control lepton masses and mixing 
angles also control the flavor pattern of the operators in Ceff- This result is interesting 
in several respects. First of all the presence of the factors can help in suppressing the 
rates of rare charged lepton decays while allowing for a relatively small and accessible scale 
A/vrp. Second, once the above expansion has been determined in a given model, it could 
be possible to establish characteristic relations among LFV processes as a consequence of 
flavor symmetries and of their pattern of symmetry breaking. Finally, if A^rp is sufficiently 
small, this opens the possibility that new particles might be produced and detected at 
the LHC, with features that could additionally confirm or reject the assumed symmetry 
pattern. All this allows, at least in principle, to realize an independent test of the flavor 
symmetry in the charged lepton sector. While the size of the scale A^vp could be relatively 
small, in our presentation we assume that the flavour scale or cutoff A is extremely large, 
possibly as large as the GUT scale. Then all low-energy effects due to the flavon dynamics 
are essentially those associated to their VEVs, which enter the effective higher dimensional 
operators through the dimensionless combination e. Virtual flavon exchanges give rise to 
other higher dimensional operators which are depleted by inverse power of A and can be 
safely neglected. A much richer variety of effects due to the flavour dynamics would be 
possible if the scale A were much smaller, close to the 100 TeV energy range, but we do 
not consider this possibility here. 



(82) 
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The effects described by C^jj are well-lcnown. In a field basis where the kinetic terms 
are canonical and the charged lepton mass matrix is diagonal the real and imaginary parts 
of the diagonal matrix elements Za are proportional to the anomalous magnetic moments 
(MDM) and to the electric dipole moments (EDM) di of charged leptons, respectively: 

Oj = 2mj^ _ „ — ReZii , di = e _ „ — ImZa . (83) 

The off-diagonal elements Zij describe the amplitudes for the radiative decays of the 
charged leptons: 

where a is the fine structure constant, Gi? is the Fermi constant and mf^ is the mass of the 
lepton /j. Finally the four-fermion operators, together with the dipole operators controlled 
by Z, describe other flavor violating processes like // — eee, r — /i/i/i, r — eee. 

An interesting example of flavor symmetry is that of minimal flavor violation (MFV) 
|241H249j whose (minimal) flavor symmetry group in the lepton sector is G/ = SU{?>)e<^ x 
SU{2t)i. Electroweak singlets e*^ and doublets I transform as (3, 1) and (1,3), respectively. 
The flavon fields or, better, their VEVs are the Yukawa couplings of the charged leptons, 
Yi = mi/v, and the adimensional coupling constants rj of the five- dimensional operator O5 
in eq. ([2]). They transform as (3,3) and (1,6), respectively. In a basis where the charged 
leptons are diagonal, we have 

= ^mf"s , rj = ^U*mf^^U^ , (85) 

where M here denotes the mass scale suppressing the operator O5. In MFV models the 
leading off-diagonal elements of Z^j are given by: 

Zij = c (Yi r]^ri)ij 

= V2c^—[AmliU.^U*^±^ml,JJ,,U*,] (86) 

where c is an overall coefficient of order one and the plus (minus) sign refers to the case 
of normal (inverted) hierarchy. We see that, due to the presence of the ratio M^/f^ the 
overall scale of these matrix elements is poorly constrained. This is due to the fact that 
MFV does not restrict the overall strength of the coupling constants r/, apart from the 
requirement that they remain in the perturbative regime. Very small or relatively large 
(but smaller than one) r] can be accommodated by adjusting the scale M. Thus, even 
after fixing A^vp close to the TeV scale, in MFV the non-observation of li — )■ Ij^ could be 
justified by choosing a small M, while a positive signal in /x — )■ 67 with a branching ratio 
in the range 1.2 x 10"^^ -t- 10"^^ could also be fitted by an appropriate M, apart from a 
small region of the ^13 angle, around ~ 0.02 where a cancellation can take place in the 
left-hand side of eq. (186|) . The dependence on the scales M and A^rp can be eliminated 
by considering ratios of branching ratios. For instance: 



2Am? 



°' ±v^ sin ^136^^ 



2 



<1 , (87) 
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where we took the TB ansatz to fix 612 and 623- We see that BR{fi — )■ 67) < BRij — )■ /i7) 
always in MFV. Moreover, for ^13 above approximately 0.07, BR{fi — > 67) < 1.2 x 10"^^ 
implies BR{t — )■ yU7) < 10~^. For ^13 below 0.07, apart possibly from a small region around 
6'i3 ^ 0.02, both the transitions /i — )■ 67 and t fi'-f might be above the sensitivity of 
the future experiments. The present limits are BR{t — )■ /i7) < 1.6 x 10~^ and BRij -> 
67) < 9.4 X 10~^. A future super B factory might improve them by about one order of 
magnitude. In the SUSY case there are two doublets in the low-energy Lagrangian and 
we should take into account the tan/3 dependence. 

A different result for the matrix Z is obtained in the model described in Sect. 4 where 
Gf = Ai X Z3 X U{1)fn- Starting from the relevant set of invariant operators, after 
the breaking of the flavor and electroweak symmetries, and after moving to a basis with 
canonical kinetic terms and diagonal mass matrix for charged leptons, we flnd |250[l251j : 

/ 0{t^e) Oi^e^) 0{t^e^) \ 
Z= 0{te^) 0{te) 0{te^) (88) 
V Oie') 0{e^) 0{e) ) 

where each matrix element is known only up to an unknown order-one dimensionless 
coefficient. There are two independent symmetry breaking parameters. The parameter 
t = {0)/A controls the charged lepton mass hierarchy and e = v^/A describes the breaking 
of A4. Notice that the uncertainty in the overall scale of the matrix elements Zij is related 
to the parameter e and is much smaller than the corresponding uncertainty in MFV. We 
can see that MDMs and EDMs arise at the first order in the parameter e. By assuming 
that the unknown coefficients have absolute values and phases of order one, from eqs. (183|) 
and (188|) we have: 

From the existing limits on MDMs and EDMs and by using eqs. fl89|) as exact equalities 
we find the results shown in table 7. 



de < 1.6 X 10 e cm 


Anp > 80 TeV 


< 2.8 X 10^^^ e cm 


Anp > 80 GeV 


5ae < 3.8 X 10-12 


Anp > 350 GeV 


Sa^ = 302 ± 88 X 10-^^ 


A^P^2.7 TeV 



Table 7: Experimental limits on lepton MDMs and EDMs and corresponding bounds on the scale Anp, 
derived from eq. (f89|) . The data on the t lepton have not been reported since they are much less 
constraining. For the anomalous magnetic moment of the muon, (5a^ stands for the deviation of the 
experimental central value from the SM expectation |2521l253j . 

Concerning the fiavor violating dipole transitions, from eq. ( I88ll we see that the dominant 
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contribution to the rate for /j — > /j7 is given by: 

\wij (90) 



BR{k Ij-i) _ AS-K^a I _ 



BR{li -> /ji^ii/j) G|,A 



4 



where are numbers of order one. As a consequence, the branching ratios of the three 
transitions /i — 67, r — /i7 and r — )■ 67 are all expected be of the same order: 

BR{ii 67) ^ BR{t /i7) ^ 5i?(r ^ 67) . (91) 

This is a distinctive feature of this class of models. Given the present experimental bound 
on BR{fi — )■ 67), eq. ( 1911) implies that r ^ fij and r — 67 have rates much below the 
present and expected future sensitivity. Moreover, from the current (future) experimental 
limit on BR{fi — )■ 67) |239y240p254] and assuming \w^e\ = 1, we derive the following bound 
on \e/A^ I- 



NP\- 



BR{n 67) < 1.2 X 10"^' (lO^^'^) 



< 1.2x10^" (1.1 X 10"^^) GeV"^ .(92) 



Taking two extreme values for the parameter |£| we find 

Anp > 20 (67) TeV [\e\ = 0.005] 

Anp > 65 (210) TeV [\£\ = 0.05] . (93) 

This model also allows for four-fermion operators that are not suppressed by any power 
of the small parameter t or e and that violate the individual lepton numbers Lj |255] . 
They are all characterized by the selection rule ALgAL^jAL^ = 2. For instance, one such 
operator is 

(/!)'([/)" = [Url^lr + Wrk + lr^el^^ + h.c] + ... . (94) 

where dots stand for additional flavor conserving contributions. These operators can 
contribute to LFV decays such as r~ — )■ fi'^e~e~, t~ — 6+/i"/i~ and their conjugate, 
whose branching ratios have upper bounds of the order of 10~^ |256] . Through a rough 
dimensional estimate we find a lower bound on the scale Aatp of the order of 15 TeV. From 
the previous considerations we see that, even invoking a cancellation in the imaginary 
part of Zee to suppress the contribution to the electron EDM, it is difficult to avoid the 
conclusion that the scale A^rp should lie considerably above the TeV range. We recall that 
if the operator in eq. (18 ip originates from one-loop diagrams via the exchange of weakly 
interacting particles of masses rriNP, then in our normalization a lower bound on A^rp of 
20 TeV corresponds to a lower bound on rriNP of about gNpANp/ [Ati) ^ 1 TeV, assuming 
(77VP similar to the SU(2) gauge coupling. 

All the previous estimates are based on an effective Lagrangian approach, with no 
explicit reference to the dynamics at the scale A^rp. If the degrees of freedom associated 
to the new physics at the scale A^p and their interactions are known, it is possible to 
directly compute the amplitudes of interest. For instance, the SUSY model of Sect. 4 can 
be completed by adding a set of soft SUSY breaking terms, which are constrained by the 
invariance under Gj = /I4 x Z3 x U{1)pm and its pattern of symmetry breaking |257H260j . 
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LFV amplitudes arise at one-loop level, via exchange of sleptons, charginos and neutralinos 
with masses of order msusY- An explicit computation of BR{li — Ij^) confirms both the 
predictions of eq. ( 191]) and the behaviour of eq. fl90|) . with A^rp = {'^t^ / g)'msusY ■ The 
coefficients Wij are typically of (9(0.1). When e is small, which also entails small tan/3 in 
our model, relatively light SUSY particles are allowed, while for e close to its upper limit, 
0.05, SUSY particle masses of several hundred GeV or close to the TeV are needed to 
satisfy the present bound on BR{fi — )■ 67), particularly if tan/3 is larger than 10. In either 
case there is only a very limited region of the parameter space where it is possible to explain 
the observed discrepancy in the muon MDM and to satisfy at the same time the current 
limit on BR{fi — )■ 67). An interesting special case is that of universal SUSY breaking 
terms, giving rise to a cancellation in the elements of Zij below the diagonal |259[I261] . 
Under these circumstances BR{li — )■ /j7) scale as rather than as e^, with the possibility 
of much lighter SUSY particles. In SUSY ^4 models also LFV 4-fermion operators are 
depleted by powers of e and the corresponding bounds on msusY are relaxed. 

In the model discussed in Sect. 8, with G j = S^x Z^x U{1)fn, the matrix Z is given 
by [262] : 

/ 0{t^v^) 0{tW) 0{tW) \ 
Z = 0{tvv') 0{tv) Oitvv'^) (95) 
V 0{vv') 0{vv'^) 0{v) ) 

Predictions for EDMs and MDMs and corresponding bounds are similar to those discussed 
above in the case of the model and summarized in Table 7. Concerning the radiative 
decays of the charged leptons we find that R^e and Rre scale as w'^/A^p, whereas Rr^ 
scales as v''^/A%p. In this case the symmetry breaking parameter v' is considerably larger 
than the parameter e of the model and this gives rise to more restrictive bounds on 
the scale of new physics A^vp. From BR{^ — )■ 67) < 1.2 x 10~^^ (10~^^) we get: 

A;vp > 90 (300) TeV [v' = 0.1] 

Anp > 130 (430) TeV [v' = 0.2] . (96) 

The model also predicts: 

BR{fx 67) ^ BR{t -> e7) > BR{t fxj) . (97) 

Summarizing, in models with discrete fiavor symmetries LFV processes are generically 
suppressed by the presence of small symmetry breaking parameters. However such a 
suppression is not completely efficient, at least in the explored models, to guarantee a 
scale of new physics close to the TeV. The best case is the one of the model, thanks to 
the very small expansion parameter e. In specific SUSY realizations of the A4 symmetry 
the present limits on the branching ratios of LFV processes still allow for a relatively light 
spectrum of superparticles, in a region of masses of interest to LHC. 

10 Leptogenesis 

The violation oi B — L implied by the see-saw mechanism suggests an interesting 
link between neutrino physics and the mechanism that produced the observed baryon 
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asymmetry in the early universe. If the baryon asymmetry originates well above the 
electroweak scale, B — L violation represents a necessary condition, since any initial B + L 
asymmetry would be erased in the subsequent evolution of the universe. According to 
leptogenesis the asymmetry is determined by the CP violating, out-of-equilibrium decays 
of the heavy RH neutrinos [263] . Through B — L non-conservation of neutrino interactions, 
the asymmetry is first generated in the leptonic number and then partly converted into 
the observed baryonic one via sphaleron interactions. Depending on whether the relevant 
decays occur at a sufficiently high temperature or not, we have an unflavored regime, where 
the leptons in the final state are indistinguishable, or a fiavored regime, where the specific 
interactions of the different leptons in the decay products cannot be neglected [26411265] . It 
is also quite remarkable that, at least in its simplest implementation, leptogenesis requires 
light neutrino masses below the eV scale |266H269] . in a range which is fully compatible 
with other experimental constraints. Unfortunately, without any additional assumptions, 
it is difficult to promote this elegant picture into a testable theory, due to the large number 
of independent parameters of the see-saw model. 

Models of lepton masses based on fiavor symmetries typically depend on a restricted 
number of parameters, thus opening the interesting possibility of relating the baryon 
asymmetry to other low-energy observables. As a general rule, to provide a realistic 
description of lepton masses and mixing angles the fiavor symmetry should always be 
broken. The breaking is described by a set of small dimensionless quantities e, which 
provide efficient expansion parameters. As we have seen in the previous sections, small 
observable quantities such as charged lepton mass ratios, 6*13, 6*23 — n/A can be expanded 
in power series of e, and the predictions are dominated by the lowest (positive) power. 

In the context of leptogenesis, given the extreme smallness of the baryon asymmetry 
[270] 

r/£^^ = (6.2±0.15) X 10"^°, (98) 

it can be convenient, at least in a certain regime, that the CP asymmetries in the RH 
neutrino decays are also suppressed by powers of e. If the baryon asymmetry is dominated 
by the decay of a single RH neutrino, we can write B 

r]B = dik (99) 

where d describes the combined effect of sphaleron conversion and dilution from photon 
production, ^ is the relevant CP asymmetry and k takes into account the wash-out effects. 
Typically we expect a dilution factor d of order 10~^ and, barring fine-tuning of the 
parameters, a wash-out factor k in the range 10^^ ^ 10^^, which favors ^ around 10^^ — 
10~^. Such CP asymmetry arises from the interference of the tree-level and the one- 
loop decay amplitudes and depends quadratically on the neutrino Yukawa couplings. In 
models like the ones discussed in Sects. 4 and 8, where the RH neutrino masses are very 
large, close to 10^^ GeV, and the corresponding neutrino Yukawa couplings are of 0(1), 
a rough estimate of the total CP asymmetry would give ^ = (9(l/(87r)), by far too large 
compared to 10^^ ^ 10~^. It is therefore interesting to analyze under which conditions 

^We will denote the CP asymmetries with ^ and we keep the letter e to indicate the generic expansion 
parameter of a spontaneously broken flavor symmetry. 
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the CP asymmetries vanish in the hmit of exact symmetry, so that the first non- vanishing 
contribution is given by some power of the symmetry breaking parameters e. If the CP 
asymmetry relevant for leptogenesis is suppressed by powers of e, this opens the very 
interesting possibihty of relating the observed baryon asymmetry rjB to other low-energy 
observable quantities [g21[2U71[2U5] such as 6^13, 6*23 - 7r/4, BR{li Ij'j). 
The total CP asymmetries in the decay of a RH neutrino are 

r - r 

6 = — — =^ (100) 

where Fj (Fj) is the decay rate of uf into leptons (antileptons). In the flavored regime the 
relevant asymmetries involve final states with a specific lepton flavor /. The flavored 
regime takes place for Mj < c 10^^ GeV where c = 1 (1 + tan^/3) in the ordinary (SUSY) 
case. The unflavored regime occurs for RH neutrino masses above that threshold. At 
one-loop we have |271] : 

OTT J/ j,- 



where 3^ is a combination of the neutrino Yukawa couplings Y = /vu 

y^J = {YY^),^ . (103) 



and the hat in eqs. f ll01lll02p denotes a basis where the mass matrix M of heavy Majorana 
neutrinos and that of charged leptons, mi, are diagonal. The functions fij and gtj depend 
on the mass ratios of the RH neutrino masses Mj. From eqs. f ll01|102p we see that both 
C,if and vanish if y is diagonal. The total asymmetries C,i vanish also if 3^ has real 
non-diagonal entries. A necessary and sufficient condition for a diagonal 3^ is: 

yM - My^ = , (104) 

where the matrices y and M are evaluated in any basis. 

If the model is invariant under the action of a flavor symmetry group Cf we have an 
interesting sufficient condition for the vanishing of the CP asymmetries. If the heavy RH 
neutrinos transform in a (three-dimensional) irreducible representation of G/, then in the 
limit of exact symmetry, where the symmetry breaking parameters e go to zero, all CP 
asymmetries vanish |272j . In this limit it is possible to show that y becomes proportional 
to the unit matrix as a consequence of a completely general group theoretical property. 
Thus, from eqs. fll01|102p we conclude that the asymmetries and vanish. Notice 
that irreducible representations of dimension larger than one are only possible if G/ is non- 
abelian. Beyond the symmetry limit, in general y gets corrections and develops complex 
off-diagonal entries at some order e^. If the spectrum of RH neutrinos is non-degenerate in 
the symmetry limit, we expect C,i = O^e'^^) and ^j/ = 0{e^). Degeneracy of RH neutrinos 
can modify this behavior through the dependence on e of the functions fij and Qij. This 
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result applies to both the models described in Sects. 4 and 8, where the RH neutrinos 
transform in the three-dimensional representions of A4 and 5*4, respectively. In the limit 
of exact flavor symmetry we find in both cases y = |yp 1 where 1 denotes the identity 
matrix. This equality holds in any basis, in particular in the mass eigenstate basis of RH 
neutrinos and we have C,i = in the symmetry limit. In both models all RH neutrino are 
very heavy, with masses well above 10^^ GeV, and the unfiavored regime applies. 

In the A4 model of Sect. 4, y acquires complex off-diagonal entries of order e ^ vt/A. 
The CP asymmetries depend only on three real parameters: two independent real 
symmetry breaking parameters Si and the lightest neutrino mass. In particular there is 
only one independent phase which is determined by the lightest neutrino mass up to an 
overall sign. We have approximately |273] 



More precisely |272p274[|275] . for normal ordering of the neutrino mass spectrum all asym- 
metries are of the same order of magnitude. For inverted ordering the two asymmetries 
^1^2 get enhanced compared to the approximate estimate of eq. (11051) by a factor 10^ 
coming from the functions fu and /21, as a result of the near degeneracy of two heavy 
RH neutrinos. To reproduce the observed baryon asymmetry, eq. (|98l) . different wash-out 
effects are required in the two cases. In the case of normal ordering the experimental value 
in eq. (^E\i is obtained when the parameter e is in its natural window, 5 x 10~^ -h 5 x 10~^, 
for a wide range of neutrino Yukawa couplings y. For inverted ordering a much larger 
wash-out suppression is needed. When e falls in the optimal range 5 x 10~'^-^5 x 10~^ this 
can be accommodated by restricting both y x sin /3 and ms in a limited range. It is quite 
remarkable that in both cases the range of the symmetry breaking parameter e suggested 
by the constraints on lepton masses and mixing angles corresponds to that required to get 
the observed baryon asymmetry through leptogenesis. 

In the 5*4 model discussed in Sect. 8, y acquires complex off-diagonal entries at the 
order v'^/v' and the CP asymmetries are expected to be of order / {v''^8tt) . Assuming 
a typical wash-out suppression of order 10~^, the observed baryon asymmetry can be 
obtained for values of {v, v') close to the range selected to fit charged lepton masses and 
mixing angles. 

Another class of models where the CP asymmetries vanish is the one of type I see- 
saw models where the Dirac and Majorana neutrino mass matrix and M as well as 
their see-saw combination rriy are form-diagonalizable. A matrix A depending on a set 
of parameters is said to be form-diagonalizable |276j if it is diagonalized by unitary 
transformations that do not depend on «j: 



where A'''[a) is diagonal and the unitary matrices Ul,r are independent from a. Examples 
of form-diagonalizable matrix are my in eqs. fll3fl4p and fll8|19p . The parameters are the 
eigenvalues mi^2,3 and the diagonalizing matrices are Ur = = Utb and Ur = Ul = 
Ubmi respectively. As we have seen in section 2, form-diagonalizable matrices naturally 
arise in the context of models with discrete flavor symmetries. It is possible to show that 




(105) 




(106) 
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if in a type I see-saw m^, M and m^, are all form-diagonalizable, then the matrix 3^ is 
diagonal and the CP asymmetries vanish |2 77(1278] . In realistic models m^, M and m,^ 
are typically form-diagonalizable only in some symmetry limit. Symmetry breaking terms 
usually spoil this property and allow for small non-vanishing CP asymmetries. 

So far we have discussed the regime of large RH masses and large neutrino Yukawa 
couplings. When the smallest RH neutrino mass is below the so-called Davidson-Ibarra 
bound |279] (4 x 10^ 2 x 10^) GeV, and we are in the regime of strong hierarchy among 
RH neutrino masses, the CP asymmetry associated to the lightest RH neutrino decay is 
too small to allow for a successful leptogenesis. To evade the Davidson-Ibarra bound we 
should depart from the strong hierarchical regime. Under certain conditions a significant 
enhancement of the CP asymmetry can be achieved even for RH neutrino mass ratios as 
small as 0.1 |280[I281] . Alternatively, we can exploit the regime of resonant leptogene- 
sis |282] , occurring when the decaying RH neutrino is quasi-degenerate in mass with some 
other RH neutrino, the mass differences being comparable with the RH neutrino decay 
width. A quasi-degeneracy of the RH neutrino spectrum is better understood and dynami- 
cally controlled in the presence of an underlying flavor symmetry. Several symmetries have 
been proposed in the literature such as Cf = SU{3)ec x SU{3)i x 0(3),yc in minimal lepton 
flavor violation pl6ti2l9ll283l[2M] or Gf = S0{3) in ref. [285]. In these two examples the 
light neutrino masses and their mixing angles are not explained but just accommodated. 
An interesting model based on a flavor symmetry group G/ = /I4 x ^3 x Z4 is that of 
ref. |286] . Like the model discussed in Sect. 4 it predicts a lepton mixing close to TB. Due 
to the presence of an additional discrete factor in the symmetry group, the RH neutrino 
spectrum is degenerate at LO, and the degeneracy is lifted by radiative corrections or 
small soft breaking terms, allowing for successful resonant leptogenesis, for a wide range 
of RH neutrino masses. 

11 Summary and conclusion 

We have reviewed the motivation, the formalism and the implications of applying non 
abelian discrete flavor groups to the theory of neutrino mixing. The data on neutrino 
mixing are by now quite precise. It is a fact that, to a precision comparable with the mea- 
surement accuracy, the TB mixing pattern is well approximated by the data (see Fig. (2)). 
If this experimental result is not a mere accident but a real indication that a dynamical 
mechanism is at work to guarantee the validity of TB mixing in the leading approximation, 
corrected by small non leading terms, then non abelian discrete flavor groups emerge as 
the main road to an understanding of this mixing pattern. Indeed the entries of the TB 
mixing matrix are clearly suggestive of "rotations" by simple, very specific angles. In fact 
the group A^, the simplest group used to explain TB mixing, is specified by the set of those 
rotations that leave a regular tetrahedron invariant. We have started by recalling some 
basic notions about finite groups and then we have concentrated on those symmetries, 
like A4 and 5*4, that are found to be the main candidates for obtaining TB mixing. We 
have discussed the general mechanism that realizes TB mixing within the framework of 
discrete flavor symmetries. The symmetry is broken down to two different subgroups in 
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the charged lepton sector and in the neutrino sector, and the mixing matrix arises from 
the mismatch between the two different residual simmetries. TB mixing requires a flavor 
symmetry group possessing appropriate residual subgroups. The breaking can be realized 
in a natural way through the specific vacuum alignments of a set of scalar flavons. We 
have described a set of models where TB mixing is indeed derived at leading order within 
this mechanism. There are many variants of such models (in particular with or without 
see-saw) with different detailed predictions for the spectrum of neutrino masses and for 
deviations from the TB values of the mixing angles. In general at NLO the different mix- 
ing angles receive corrections of the same order of magnitude, which are constrained to 
be small due to the experimental results which are very close to the TB values. Indeed 
the small experimental error on 612, which nicely agrees with the value predicted by TB 
mixing, suggests that the NLO corrections should be of order of few percent, at most. 
Additional symmetries are needed, typically of the U{1)fn or type, in order to repro- 
duce the mass hierarchy of charged leptons. In the neutrino sector there is no reason for 
the mass eigenvalues not to be of the same order in absolute value. Thus the smallness of 
the ratio ^/r ~ 0.2, where r is defined in eq. ([1]), is accidental in most of these models. 
Both normal and inverse hierarchy spectra can be realized. The phenomenology of the 
models was summarized. We have also discussed the implications of models based on dis- 
crete flavor groups for lepton flavor violation and for leptogenesis. Lepton flavor violating 
processes, the muon g-2 and the EDM's of leptons impose strong constraints on every new 
physics model. This is also true for the models considered here. But the specific suppres- 
sion factors and selection rules induced by the finite flavor symmetry group, in particular 
by A4, may help to improve the consistency of the model even in the presence of new 
physics at the TeV scale. The observed baryon asymmetry in the Universe, explained in 
terms of leptogenesis from the decay of heavy Majorana neutrinos, is found to be com- 
patible with models based on discrete groups. Neutrino Yukawa couplings of order one 
and RH neutrino masses of order 10^^ -^ 10^^ GeV would typically lead to CP asymmetries 
too large to reproduce the observed baryon asymmetry. However, as a consequence of a 
general group theoretical property, in all models where the three RH neutrinos transform 
in a single irreducible representation of the flavor group, the unflavored CP asymmetries 
vanish in the limit of exact symmetry and small values can be generated through NLO 
corrections. 

An obvious question is whether some additional indication for discrete flavor groups can 
be obtained by considering the extension of the models to the quark sector, perhaps in a 
Grand Unified context. The answer appears to be that, while the quark masses and mixings 
can indeed be reproduced in models where TB mixing is realized in the leptonic sector 
through the action of discrete groups, there are no specific additional hints in favour of 
discrete groups that come from the quark sector. Examples of Grand Unified descriptions 
of all fermion masses and mixings with TB mixing for neutrinos have been produced and 
have been discussed in this review. For quarks, only the third generation masses are present 
at leading order in these models. The other entries of the mass matrices are small due to 
additional symmetries or other dynamical reasons (for example, suppression factors from 
extra dimensions), and the small mass ratios and the small mixing angles are generated 
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by these corrective effects and are not due to the discrete group. As a consequence, the 
action of the discrete ffavor group is only clearly manifest among the comparable neutrino 
sector masses, in the basis where charged leptons are diagonal. 

Different forms of neutrino mixing other than TB mixing are also amenable to a de- 
scription in terms of discrete groups. In alternative to TB mixing, in sect. 8 we have 
discussed the possibility that actually a more appropriate starting point, could be BM 
mixing, corrected by large terms of 0(Ac), with Ac being the Cabibbo angle ("weak com- 
plementarity"), arising from the diagonalization of charged leptons. By suitably modifying 
the construction in terms of discrete groups adopted in the case of TB mixing, we have 
identified the group 5*4 as a good candidate to also provide, in a different presentation, the 
basis for naturally obtaining BM mixing in first approximation. In the model described 
the NLO terms are such that the dominant corrections only affect 6*12 and 6*13 (which re- 
ceive 0{Xc) shifts), while 623 receives smaller corrections. A value of 6*13 near the present 
bound would support this possibility. 

In the near future the improved experimental precision on neutrino mixing angles, in 
particular on ^13, could make the case for TB mixing stronger and then, as a consequence, 
also the case for discrete flavor groups would be strenghtened. Further important input 
could come from the LHC In fact, new physics at the weak scale could have important 
feedback on the physics of neutrino masses and mixing. 
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